THE FLORIDA STATE UNIVERSITY 
COLLEGE OF ARTS AND SCIENCES 



NEXT-TO-LEADING-ORDER CORRECTIONS TO WEAK BOSON 
PRODUCTION WITH A MASSIVE QUARK JET PAIR AT HADRON 

COLLIDERS 



By 

FERNANDO FEBRES CORDERO 



A Dissertation submitted to the 
Department of Physics 
in partial fulfihment of the 
requirements for the degree of 
Doctor of Philosophy 



Degree Awarded: 
Fall Semester, 2007 



a Kiko y Tefina... 



ii 



TABLE OF CONTENTS 



3. Numerical Results 

3.1 The Setup 

3.2 Wbb Production at the Tevatron 

3.3 Zbb Production at the Tevatron . 



4. CONCLUSION 
4.1 Outlook . 



APPENDICES 

A. Standard Model of Particle Physics 



List of Tables 

List of Figures 

Abstract Q 

1. INTRODUCTION [H 

2. NLO Calculation for W/Z bb Production at Hadron Colliders 

2.1 Factorization Theorem for Hadronic Cross Sections 

2.2 Hadronic Cross Section for W/Z bb Production: General Structure at NLO 

2.2.1 Renormalization Scale Dependence of the NLO Cross Section . . . 

2.2.2 Calculating the Virtual Cross Section a™^ 

2.2.3 Calculating the Real Cross Section a'^^^^ 

2.3 NLO QCD Corrections to Wbb Production at Hadron Colliders 

2.3.1 Virtual Corrections to qq' —>■ Wbb 

2.3.2 Real Corrections to Wbb Production 

2.3.3 Total Cross Section for pp{pp) Wbb 

2.4 NLO QCD Corrections to Zbb Production at Hadron Colliders 

2.4.1 Tree Level Cross Section for Zbb Hadronic Production 

2.4.2 Virtual Corrections to qq — > Zbb 

2.4.3 Virtual Corrections to gg Zbb 

2.4.4 Real Corrections to Zbb Production 

2.4.5 Total Cross Section of pp{pp) Zbb 



iii 



B. Tree Level Amplitudes for W/Z bb Production 

C. Scalar Integrals for W/Z bb Production 

D. Reducing Tensor Feynman Integrals 

E. Using Quadruple Unitarity Cuts to Check Coefficients of Scalar Box Integrals 

F. Phase Space Integrals for the Emission of a Soft Gluon in the two Cutoff PSS 
Method 

REFERENCES 



iv 



LIST OF TABLES 



Diagram content of color amplitudes and whose color factor, when 
interfered with Aq'' or Aq""^ respectively, is Ci 

Diagram of content color amplitudes Af' and A^"''' whose color factor, when 
interfered with Aq'^ or Aq""^ respectively, is C2 

Diagram content of color amplitude Af' whose color factor, when interfered 
with Aq^ or Aq"''' respectively, is C3 

Diagram content of color amplitude A^""^ whose color factor, when interfered 
with Aq'' or Aq'''' respectively, is NC2 

Diagram content of color amplitude A^'' whose color factor, when interfered 
with Af or respectively, is (A^ - 2/N)C2 

LO and NLO total Wbb cross sections at the Tevatron for massive and massless 
bottom quarks, using fir = fJ'f = Mw + 2mf). The numbers in square brackets 
are the ratios of the NLO and LO cross sections, the so called i^-factors. 
Statistical errors of the MC integration amount to about 0.1% 

LO and NLO total Zbb cross sections at the Tevatron for massive and massless 
bottom quarks, using fi^ = fJ^f = Mz + 2mb. The numbers in square brackets 
are the ratios of the NLO and LO cross sections, the so called i^'-factors. 
Statistical errors of the MC integration amount to about 0.1% 



V 



LIST OF FIGURES 



1.1 Example of processes to which Vbb (V = W or Z) is a main background. . . 

2.1 Quadruple cut check of the calculation of a box diagram involving a top-quark 
loop [l5 

2.2 Tree level Feynman diagrams for qq' Wbb 



fl 2,') fl 2l 

2.3 Gluon {S[ ' ) and quark (5*2 ' ) 0{as) self-energy corrections contributing 
to the qq' — >• Wbb subprocess at NLO. The shaded blobs denote standard 



2.4 0{as) vertex corrections contributing to the subprocess qq' Wbb at NLO. 
The shaded blobs denote standard one-loop QCD corrections to the qqg, bbg 



2.5 0{as) box diagram corrections contributing to the qq' Wbb process at 

(1 2) 

NLO. The term Crossed refers to the box diagrams, B^^' obtained from the 

>(1,2) 



2.6 0{as) pentagon diagram corrections contributing to the qq' Wbb process 



2.7 0{as) real corrections: examples of initial and final real gluon emission and 



one-loop QCD corrections to the gluon and quark propagators, respectively. . 120 



and qq'W vertices, respectively |21 



i?3 ' boxes by flipping the 6-quark fermion line |22 



at NLO 22 



q{q)g initiated subprocess |26 



2.8 Dependence of the total Wbb NLO QCD cross section on the 6s PSS parameter. |34 



2.9 Dependence of the total Wbb NLO QCD cross section on the 6c PSS parameter. |35 



2.10 Tree level Feynman diagrams for qq — ^ Zbb, with Z emitted from final fermion 



line 36 



2.11 Tree level Feynman diagrams for gg — >■ Zbb. The circled crosses correspond to 



all possible insertions of the Z boson, each one representing a different diagram. |36 



vi 



(1 2") (1 2') 

2.12 Gluon {S\X') and 6-quark {S'^^f') C>(a,) self energy corrections contributing to 
the qq — * Zbb subprocess at NLO, when the Z boson is emitted from the final 
fermion line. The shaded blobs denote standard one-loop QCD corrections to 
the gluon and quark propagators respectively. 

2.13 0{as) vertex corrections contributing to the qq Zbb subprocess at NLO 
when the Z boson is emitted from the final fermion line or from a closed 
fermion line. The shaded blobs denote standard one-loop QCD corrections to 
the qqg (V^/.f'^^), bbg (V^^''^^) and qqZ {V^^^''^^) vertices respectively. V^*'.]-'^'* are 
b- and t-fermion loop vertices which are UV and IR finite (contributions of 
quarks from first and second family vanish) 

2.14 0{as) box diagram corrections contributing to the qq Zbb subprocess at 
NLO, when the Z boson emitted from the final fermion lines {b or b) 

2.15 0{as) pentagon diagram corrections contributing to the qq — > Zbb subprocess 
at NLO, when the Z boson emitted from the final fermion lines {b or b). . . . 

2.16 Gluon {S[]f^) and quark {S[];!^\ Si]f\ S^];^ly and S^^^l^^^f^) self 
energy corrections contributing to the gg Zbb subprocess at NLO. The 
circled crosses correspond to all possible insertions of the Z boson, each one 
representing a different diagram 

2.17 0{as) vertex corrections contributing to the gg Zbb subprocess at NLO. 
The shaded blobs denote standard one-loop QCD corrections to the ggg 
(f/;/^ bbg and V^}ff^) and bbZ (V^^'^ ^ and V,^};^f^) vertices. 
The circled crosses correspond to all possible insertions of the Z boson, each 
one representing a different diagram 

2.18 0{as) box diagram corrections contributing to the gg — > Zbb subprocess at 
NLO. The circled crosses correspond to all possible insertions of the Z boson, 
each one representing a different diagram 

2.19 0{as) pentagon diagram corrections contributing to the gg — > Zbb subprocess 
at NLO. The circled crosses correspond to all possible insertions of the Z 
boson, each one representing a different diagram 

2.20 Examples of 0{as) real corrections to qq Zbb production, with Z emitted 
from the final fermion line 

2.21 Examples of 0{as) real corrections to Zbb production. The circled crosses 
denote all possible insertions of a Z weak boson, each insertion corresponding 
to a different diagram 

2.22 Dependence of the total Zbb NLO QCD cross section on the 6s PSS parameter 

2.23 Dependence of the total Zbb NLO QCD cross section on the 6c PSS parameter 



vii 



3.1 Dependence of the LO (black solid band), NLO inclusive (blue dashed 
band), and NLO exclusive (red dotted band) Wbb total cross sections on the 
renormalization/factorization scales, including full bottom-quark mass effects. 
The bands are obtained by independently varying both fir and fif between 
/io/2 and 4/io (with fiQ = nib + Mw /'^) 

3.2 Dependence of the LO and NLO inclusive and exclusive pp — > Wbb total cross 
section on the renormalization/factorization scale, when /^^ = /xj = yU. The 
left hand side plot compares both LO and NLO total cross sections for the 
case in which the bottom quark is treated as massless (MCFM) or massive 
(our calculation). The right hand side plot shows separately, for the massive 
case only, the scale dependence of the qq' and qg + qg contributions, as well 
as their sum 

3.3 Dependence on the renormalization/factorization scale of the rescaled differ- 
ence between our NLO calculation (with ^ 0) of the total Wbb cross 
section and the corresponding result computed using MCFM (with nib = 0) 
for the inclusive and exclusive cases (with fij. = fif = f^o) respectively. The 
error bars indicate the statistical uncertainty of the Monte Carlo integration. 

3.4 The distribution da{pp — > Wbb)/dnibb in LO and NLO QCD. The right hand 
side plot shows the ratio of the LO and NLO distributions 

3.5 The inclusive and exclusive distributions da{pp Whb) / dntbi derived from 
our calculation (with mb 7^ 0) and from MCFM (with nib = 0). The right 
hand side plot shows the ratio of the two distributions, da{nib 7^ 0)/da{nib = 0). 

3.6 The LO distribution da{pp — > Wbb)/dnibb derived from our calculation (with 
nib 7^ 0) and from MCFM (with nib = 0). The right hand side plot shows the 
ratio of the two distributions, da{mb 7^ Q)/da{mb = 0) 

3.7 The mbb distribution of the rescaled difference between our NLO calculation 
(with mb 7^ 0) and MCFM (with ntb = 0) for the inclusive (upper plot) and 
exclusive (lower plot) pp Wbb production 

3.8 Dependence of the LO (black solid band), NLO inclusive (blue dashed 
band), and NLO exclusive (red dotted band) Zbb total cross sections on the 
renormalization/factorization scales, including full bottom-quark mass effects. 
The bands are obtained by independently varying both fi^- and fif between 
/io/2 and 4/io (with fiQ = nib + Mz/I) 

3.9 Dependence of the LO and NLO inclusive and exclusive pp Zbb total 
cross section on the renormalization/factorization scale, when /t^ = /i/ = /i. 
The LHS plots compare both LO and NLO total cross sections for the case 
in which the bottom quark is treated as massless (MCFM) or massive (our 
calculation). The RHS plots show separately, for the massive case only, the 
scale dependence of the qq, gg and qg + qg contributions, as well as their sum. 



viii 



3.10 Dependence on the renormalization/factorization scale of the rescaled differ- 
ence between our NLO calculation (with 7^ 0) of the total Zbb cross section 
and the corresponding result computed using MCFM (with mf, = 0) for the 
inclusive and exclusive cases (with fir = fif = fio) respectively. The error bars 
indicate the statistical uncertainty of the Monte Carlo integration 

3.11 The distribution da{pp Zbb)/dmfji in LO and NLO QCD. The right hand 
side plot shows the ratio of the LO and NLO distributions 

3.12 The inclusive and exclusive distributions da{pp Zbb)/dm^^i derived from 
our calculation (with mi, 7^ 0) and from MCFM (with = 0). The right 
hand side plot shows the ratio of the two distributions, da{mb 7^ 0)/do-{mh = 0) 

3.13 The LO distribution da{pp Zhb) / dnii^i derived from our calculation (with 

7^ 0) and from MCFM (with = 0). The right hand side plot shows the 
ratio of the two distributions, da{mi, 7^ 0)/da{'mi, = 0) 

3.14 The rrif^i distribution of the rescaled difference between our NLO calculation 
(with nib 7^ 0) and MCFM (with rrib = 0) for the inclusive (upper plot) and 
exclusive (lower plot) pp Zbb production 

C. l Topology of one-loop Feynman integrals. We denote by t the loop momentum, 

by {qi} the set of incoming momenta and by {rrii} the set of masses in the 
internal propagators 

D. l Topology of one- loop Feynman Integrals. We denote by t the loop momentum, 

{qi} the set of incoming momenta and {rrii} the set of masses in the internal 
propagators 

E. l Topology of the example presented. It corresponds to two Feynman Diagrams 

given by the two possible orientations of the fermion line. The loop in the 
fermion is a top quark (with mass rrit) 



ix 



ABSTRACT 



We present the calculation of Next-to-Leading-Order Quantum Chromo Dynamics cor- 
rections for the production oi &W oi Z weak boson associated with a bottom anti-bottom 
quark pair at hadron colliders (pp, pp W/Z bb), including the effects of a non-zero bottom- 
quark mass. We find a considerable reduction of the renormalization and factorization scale 
dependence of our results with respect to Leading-Order calculations. In particular, we study 
the impact of the corrections on the total cross section and invariant mass distributions of the 
bottom anti-bottom quark pair at the Fermilab Tevatron pp collider. We perform a detailed 
comparison with a calculation that considers massless bottom quarks and find significant 
deviations in regions of phase space with small invariant mass of the bottom anti-bottom 
quark pair. 

Our results will be relevant to ongoing and future searches at hadron colliders, as the 
W/Z bb production mode is the main background to important signals, such as light Standard 
Model Higgs boson production or single top-quark production. 
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CHAPTER 1 



INTRODUCTION 



The Standard Model (SM) of particle physics is today the best mathematical framework 
to understand the dynamics of all elementary particles that have been observed in high- 
energy collisions. This quantum field theory has been amazingly successful in describing 
and predicting to an unprecedented level of precision observations made in all high energy 
particle accelerator experiments. 

Extraordinary phenomenological efforts in the past few decades, both from the experi- 
mental side producing impressive amounts of precision data, and from the theoretical side 
producing powerful techniques to handle challenging calculations, show the SM to be a robust 
model of fundamental interactions. So much so that, if we consider only collider data, the 
SM fits all observables, with only a few showing statistical deviations of up to 2 — 3 cr q 

At the same time, however, we know that the SM is an incomplete theory, as it does not 
include gravitational interactions. Moreover, from cosmological data we know that the SM 
falls short of explaining the origin of Dark Matter and Dark Energy, and does not predict 
as large an asymmetry between matter and antimatter as observed in the universe. Another 
puzzle is the mechanism that breaks the electroweak SU{2)w x f/(l)y symmetry. In the 
SM this breaking occurs as a doublet of complex scalar fields acquires a non-zero vacuum 
expectation value (VEV), thus producing what is called spontaneous symmetry breaking. 
This mechanism gives rise to effective mass terms for weak force carriers, quarks, and leptons 
and leaves a physical scalar particle, the Higgs particle, which has so far eluded observatioEC- 
Among others, these are the reasons why nowadays there are many models that posit physics 
beyond the SM, which often embed different mechanisms of electroweak symmetry breaking 

^See for example the LEPEWWG website at http://lcpewwg.web.cern.ch/LEPEWWG and the reviews 
on SM related topics of the PDG [1] 

^ A more detailed description of this mechanism and a brief introduction to the SM is given in Appendix El 
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(EWSB). 

The current and future hadron colliders, i.e. the Tevatron, a pp collider currently taking 
data at 1.96 TeV center of mass energy at Fermilab, near Chicago, and the Large Hadron 
Collider (LHC), a pp collider with 14 TeV center of mass energy that will start in 2008 at 
CERN, Geneva, have as a main goal the elucidation of the mechanism of EWSB as well as 
the exploration of the energy spectrum beyond the weak scale, where physics beyond the 
SM (BSM) is expected. The processes studied in this thesis, i.e. the hadronic production of 
a weak force carrier with a bottom-antibottom quark pair {pp,pp — > W/Z bb), play a crucial 
role in some of the current studies of EWSB and BSM. They represent an important QCD 
background in the searches for a light SM-like Higgs boson (H) and for single top-quark 
production. 

The Necessity of Higher Order Corrections in perturbative QCD 

In order to improve our understanding of the behavior of fundamental particles at high 
energies, theorists are faced with the necessity to calculate signal and, often, background 
processes with high precision. This last task becomes essential when the signal to background 
ratio is small and the background cannot be easily extracted from data. Typically, this is 
the case for processes that involve a large number of kinematic variables and that have broad 
kinematic distributions, as often arises when final states consist of several jets and/or missing 
energy. 

At hadron colliders, QCD effects are particularly important and must be taken into 
account to obtain precise theoretical predictions. Since at high energies QCD is a pertur- 
bative quantum field theory (pQFT), QCD effects at collider energies can be calculated 
order by order in the strong coupling constant. The lowest order at which a process can be 
calculated, the Leading Order (LO), typically has a large theoretical uncertainty associated 
with it. This is mainly due to the opening of new production channels at higher orders of the 
perturbative series and to the large dependence of LO calculations on renormalization and 
factorization scales, in certain renormalization prescriptions. Adding the first order QCD 
corrections, that is, Next-to-Leading-Order (NLO) corrections, usually improves the stability 
of theoretical predictions considerably and tests the behavior of the perturbative expansion. 
Occasionally, when the NLO corrections are unusually large, the reliability of the predictions 
can be improved by computing Next-to-Next-to-Leading Order (NNLO) QCD corrections. 
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In all known examples this is enough to reduce the theoretical uncertainty to an acceptable 
level. 

Today the standard for hard scattering cross section calculations is NLO, and since the 
early nineties a wide set of processes have been studied at this level in perturbation theory 
(for an up-to-date review of some of them and a look at state of the art techniques that have 
been developed see Ref. [2]). 

However, while several programs exist that allow automated calculations of partonic 
differential cross sections at LO (e.g. Madgraph [31 HE], CompHEP [H], AMEGIC++ [7]), 
there are as yet no algorithms able to deal with all processes at NLO in a completely 
automatic way. There is a bottleneck that occurs with the calculation of virtual one-loop 
diagrams with many external partons or "legs" . Basically, following a traditional Feynman 
diagram approach, the complexity of the analytic expressions grows exponentially with the 
number of legs and the number of massive internal/external particles. This is mainly due 
to the number of Feynman diagrams, the increased number of kinematic variables and the 
increased complexity of the tensor integrals appearing in each Feynman diagram. 

It has been observed that certain amplitudes (for example the so called Maximally 
Helicity Violating (MHV) amplitudes) show a surprising analytic simplicity, hidden by the 
cumbersome intermediate steps of standard calculations. A good part of the progress in the 
field is due to a better understanding of such amplitudes (the literature on MHV amplitudes 
is now extensive, but for a brief review see Refs. [H [9l [T0|, [H]). Many new techniques have 
appeared that exploit general properties of gauge field theories such as gauge invariance, 
factorization, unitarity and the existence of representations in terms of Feynman integrals [lOl 
[T2| [21 [13] . An example of their success is the recently completed one-loop calculation of the 
set of all helicity amplitudes with six external gluons [Ill[l5l[l6l[I71[T8l[l9l[20l[2Tl[22]. We 
have used generalized unitarity, specifically quadruple cuts [23], as a means to make non- 
trivial cross checks of coefficients of scalar box integrals of sets of Feynman diagrams. This 
represents the first direct application of this techniques to a phenomenologically relevant 
computation including massive internal and external particles. 

In this dissertation we present the calculation of NLO QCD corrections to the production 
of a.W OT Z weak gauge boson in association with a bottom-antibottom quark pair at hadron 
colliders {pp,pp — W/Z bb), including full bottom-quark mass effects. The main difficulty we 
encounter is the calculation of virtual one-loop diagrams with up to five legs that include the 
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(a) Associated VH production 



(b) Single top-quark production 



Figure 1.1: Example of processes to which Vbb {V = W ot Z) is a main background. 

full effects of massive bottom-quarks. The latter increases the complexity of the calculation 
due to the addition of an extra kinematic invariant. We follow a traditional Feynman diagram 
approach to produce fully analytical expressions for one-loop amplitudes, which will allow 
non-trivial cross checks with on-shell recursion techniques, beyond the box coefficients checks 
presented in this work. 

W/Z bb Production at Hadron Colliders 

The associated production of a W/Z boson with a bb pair is by itself an interesting 
signature for hadron colliders, since it can be precisely studied experimentally. Sophisticated 
techniques exist to detect weak bosons at hadron colliders, especially when they decay into 
leptons. Even in the messy environment of a high energy hadron collider, Drell-Yan processes 
{W, Z or 7 production) are often used as tools for detector calibration and luminosity 
measurements. Furthermore, 6-quarks are very useful tools too, as the mesons in which 
they fragment have a life-time long enough to allow tagging. Nowadays the efficiency of this 
tagging is close to 50% for a large range of transverse momenta of the 6-jets. These sorts of 
studies, in the context of W/Z bb production, allow non-trivial cross checks of experimental 
techniques, and provide further constraints on the SM, as long as the theoretical uncertainty 
on such signals is equal or smaller than the experimental precision. 

The interest in W/Z bb production at hadron colliders is increased by the characteristic 
features of a light SM Higgs boson H. One of the strengths of the SM Higgs sector is that it 
is highly predictive as only one of its parameters is not yet constrained by direct observation. 
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the Higgs boson mass rriH. Electroweak precision measurements hint at the existence of a 
hght Higgs boson, with mass below 144 GeV at 95% confidence level [f. For such a Higgs 
boson one of the main production channels is the W/Z H associated production, with the 



Higgs boson decaying most of the time into a hh pair, as depicted in Figure 1.1(a) , for which 
W/Z hh represents a major background. 

The production of a Higgs boson in association with an electroweak gauge boson, 
pp HV {V = Z, W) with H hh, is indeed the most sensitive production channel of a SM 
Higgs boson at the Tevatron for a Higgs boson lighter than about 140 GeV [2 ^ [25 | [26| [27 1 [28] . 
The Tevatron with an integrated luminosity of 6 fb~^ will be able to exclude a Higgs boson 
with 115 GeV < Mh < 180 GeV at 95% confidence level [29], and will provide crucial 
guidance for the search strategy at the LHC. 

On the other hand, due to the huge hadronic activity, the experiments at the CERN LHC 
will probably first look for a light SM-like Higgs boson in the if 77 decay channel. In spite 
of the small branching ratio, ^ 77 produces a clear peak in the invariant mass distribution 
of the pair of photons. In order to fully identify a potential Higgs boson candidate, however, 
the LHC experiments will have to measure H ^ hh and this will have to be done when the 
Higgs boson is produced via associated production, either W/ZH or tiH. 

Finally, Whb is among the most relevant irreducible background processes for single-top 



quark production [301 (SH [22], as illustrated in Figure 1.1(b), both at the Tevatron and at 
the LHC. This is particularly relevant to the Tevatron, where, via single-top production, the 
Wht vertex is being measured for the first time [551 [M] . 

The cross section for pp HV has been calculated including up to NNLO QCD 
corrections [351 ESI [37] and 0(a) electroweak corrections [38], while single-top production 
has been calculated at NLO in QCD [391 [iQl [HI [121 [l3l [H [13 [ig [S], and at one-loop 
of electroweak (SM and MSSM) corrections [H]. Thus, to fully exploit the Tevatron's and 
LHC's potentials to detect the SM Higgs boson and to impose limits on its mass, as well as 
to test the third generation quark coupling to the W boson, it is crucial that the dominant 
background processes are also precisely calculated. 

In the present experimental analyseqj, the effects of NLO QCD corrections on the total 
cross section and the dijet invariant mass distribution of the W/Z hh background process 



■^For an update see the LEPEWWG website at http://lepcwwg.wcb.ccrii.ch/LEPEWWG 
^For updated results. sec the CDF and DO websites at http://www- 
cdf.fnal.gov/physics/exotic/exotic.html and http://www-dO.fnal.gov/Run2Physics/WWW/results/higgs.htni. 
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have been taken into account by using the MCFM package [H]. In MCFM, the NLO 
QCD predictions of both total and differential cross sections for the pp{p) — * W/Z bb 
production processes have been calculated in the zero bottom-quark mass (m^ = 0) 
approximation [SUl EH \S2\, using the analytical results of EH- From a study of the 
Leading Order (LO) cross section, finite bottom-quark mass effects are expected to affect 
both the total and differential W/Z bb cross sections mostly in the region of small 66-pair 
invariant masses [S2]. Given the variety of experimental analyses involved both in the search 
for HW/Z associated production and single-top production, it is important to assess precisely 
the impact of a finite bottom-quark mass over the entire kinematical reach of the process, 
including complete NLO QCD corrections. 

Using the MCFM package [H], we compared our results with the corresponding results 
obtained in the mb = limit. Numerical results are presented for the total cross section and 
the invariant mass distribution of the bb jet pair {m^i), at the Tevatron pp collider, including 
kinematic cuts and a jet-finding algorithm. In particular, we apply the /cr jet algorithm and 
require two tagged 6-jets in the final state. 

We have found that the NLO QCD corrections reduce considerably the theoretical 
uncertainty in the total hadronic cross section for W/Z bb production. We have found that 
NLO corrections are significant. The shape of the m^i invariant mass distribution is changed 
by the NLO QCD corrections; that is, the effect does not amount to a simple NLO/LO 
rescaling factor {K-i&ctoi). Finally, we have found that mass effects affect the total cross 
section by about 8% to 10% at NLO. The influence is greatest in the small m^^ invariant 
mass region. 

This dissertation is organized as follows. In Chapter [2] we present all the details of the 
full NLO calculation, including a thorough description of the techniques used to obtain both 
virtual and real corrections, first for Wbb production and then for Zbb production. We 
have included in the Appendices RllFl a few complementary reviews and several collections of 
technical details, omitted for brevity and aesthetic reasons from the body of this dissertation. 
In Chapter Owe present the results of our calculation. We show the LO and NLO dependence 
on renormalization and factorization scales, we study the impact of the corrections on the 
total cross sections and invariant mass distributions of the bottom quark-antiquark pair. We 
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show the rrih effects by comparing to a calculation that considers massless bottom-quarks. 
In Chapter H] we conclude with a summary of the main results and discuss future studies 
of the W/Z bb production mode, as well as the natural generalization of our calculation to 
other processes. 
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CHAPTER 2 



NLO Calculation for W/Z bb Production at Hadron 

Colliders 

In this Chapter we present the details of the NLO QCD calculation of Wbb and Zbb hadronic 
production including full 6-quark mass effects. 

The NLO QCD calculation of the W/Z bb production cross section for massless b quarks 
has been available in the literature for quite some time [501 EB 122] • It was performed by 
using the analytic expression of the scattering amplitudes for a weak gauge boson to four 
(massless) partons [531 El], and simulating 6-quark mass effects by imposing the kinematic 
conditions: 

{pb + pif>AQ\ pI>Q, pI>Q, (2.1) 

where Q is a scale of the order of the 6-quark mass, Pf,(i^ is the momentum of the b{b) quark 
and p'^g^^ represents the 6(6)-quark transverse momentum. 

We have improved on the massless calculation by considering a fully massive b quark both 
at the level of the scattering amplitude and in the integration over the final state phase space. 
We keep the weak bosons as on-shell particles, though the extension to include their leptonic 
decays does not present in principle any special complications. The rest of this Chapter is 
organized as follows. In Section [?!T] we present the main theoretical framework that allows the 
calculation of hadronic cross sections. We introduce briefly the parton model and the QCD 
factorization theorems. In Section !?^ we present the main characteristics of the cross sections 
for W/ Z bb hadronic production. Sections 12.31 and 12.41 present the core of the calculation for 
Wbb and Zbb production at hadron colliders respectively. They are organized into subsections 
which discuss the LO and the NLO virtual and real corrections of the calculation. We show 
explicitly the cancellation of UV singularities, by renormalization, and of IR singularities, by 
matching virtual and real corrections and by consistently absorbing all long-distant physics 
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into the renormalized parton distribution functions (PDFs) (see Section [TT]) . 

Due to the complexity of this calculation, all results have been cross checked using at least 
two independent sets of codes. The analytic calculation of the scattering amplitudes has been 
implemented using, at different stages, FORM [SS], TRACER [^S], Maple and Mathematica. 
Final numerical results have been obtained with codes built in C and FORTRAN, and we 
have used the FF package [^Z] and Madgraph [31I11[S] to cross check pieces of our code. Some 
of the figures in this Chapter have been produced using AXODRAW [SB], FeynDiagram 
and Grace [5U] . 

2.1 Factorization Theorem for Hadronic Cross 

Sections 

The parton model of hadron structure [HH E2] paved the way to the full formulation of QCD, 
as it created a framework to connect data from high energy hadronic collisions with the 
quark model according to which quarks are the elementary constituents of hadrons. The 
parton model basically states that when hadrons interact via a high momentum transfer, 
they appear as built of point-like, quasi-free constituents, the partons which we now know 
as the quarks and gluons of QCD. To each parton z in a given hadron A is associated a 
parton density function (PDF), J-'f{x), which describes the probability of finding parton i 
in hadron A with a fraction x of the total momentum of the hadron (for more details on the 
development of the model see for example Ref. [621 [63l [M]). 

Despite its success, the parton model remained without a firm theoretical basis until 
the 80's when a series of QCD factorization theorems rigorously proved that most hadronic 
observables can be calculated in QCD by disentangling the non-perturbative properties of 
hadrons and the asymptotically free behavior of partons into well defined building blockj^ 
(see for example Ref. [65l [68l [69]). As a result, a given observable Oa{Q) involving an 
initial hadronic state A and a momentum transfer Q considerably larger than Aqcd can be 
decomposed as: 

Oa{Q) = Yl fAiQ, /^/) ® 0,{Q, ^/) , (2.2) 

i 

^Bascd on general physics principles these theorems have been extended to a wide variety of processes, 
although full proofs only exist for a few processes [S5] . Recently it has been shown [BSl [ST] that certain 
processes in colliders with unpolarized beams can break factorization at Next-to-Next-to-Next-to-Leading 
Order (NNNLO), an order still far from being tested. 
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where the sum is over a full set of partons i, /i/ is the so called factorization scale and 
the operation ® is an integral convolution. The Oi{Q,fif) are point-like operators that 
describe the short-distance parton interactions which, as the momentum transfer is large, 
can be calculated perturbatively. The functions f\{Q,jjif) contain all the long-distance 
physics information related to hadron A. They are universal (i.e. process independent) 
non-perturbative objects whose Q-evolution can be determined perturbatively starting from 
a data-driven non-perturbative core. 

The factorization scale yU/ is introduced to define a boundary between the perturbative 
and non-perturbative regimes. Although the LHS of Eq. (12.21) is in principle independent of 
/i/, this is only true when considering all orders in the perturbative expansion. In practice, 
the perturbative series is always truncated and that leaves a spurious /i/ dependence in the 
calculated Oa{Q)i which can be argued that asymptotically is of the order of the next order 
in the perturbative expansion. That is why when a fixed-order calculation is performed in 
perturbative QCD, the spurious dependence on the factorization scale (and similarly on the 
renormalization scale) is used as an indicator of the theoretical uncertainty associated with 
the calculation. For the reader interested in expanding on this subject, see for example the 
review papers in Ref. [Ml [69l [70] and references therein. 

2.2 Hadronic Cross Section for W/Z hh Production: 

General Structure at NLO 

Enforcing the factorization properties of QCD cross sections, we can write the NLO QCD 
total or differential cross section for pp{pp) — > W/ Z bb as: 



where J^p^' are the PDFs for parton i in a proton (antiproton). The sum runs over all relevant 
subprocesses contributing to the hadronic differential cross section initiated by partons i and 
j. The partonic cross section for the subprocess ij — W/Z bb is denoted by Cij. The hadronic 
process pp{pp) Wbb receives contributions from the initial state qq' at LO (where q and q' 
represent quarks of up-type and down-type respectively), and from qq', qg and qg at NLO, 
while pp{pp) — s> Zhb receives contributions from qq and gg at LO and from qq, gg, qg and qg 




^^^{xi, fi)J^f^\x2, fi)aij{xi, X2, fJ.) + {xi ^ X2) , (2.3) 
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at NLO. The scale /i corresponds to both the factorization scale (/ij) and renormalization 
scale (/ir)- The factor in front of the integral is a symmetry factor that accounts for the 
presence of identical particles in the initial state of a given subprocess {6ij is the Kronecker 
delta). The partonic center-of-mass energy squared, s, is given in terms of the hadronic 
center of mass energy square, sh, by s = XiX2Sh- 

We then write the NLO partonic cross section as follows: 

= alf{x,,X2,fi) + Saf''{x^,X2,^i) , (2.4) 

where as{f^) is the strong coupling constant evaluated at the scale fi, a}f'{xi,X2, fi) is the 
0(0^) LO partonic cross section and 6afj^'~'{xi,X2, fi) contains the 0{as) corrections to 
the partonic LO cross section. The 6afj^'~'{xi,X2, ^) corrections can be decomposed in the 
following way: 

Saf'' = j d {PS;) J]| ^ W/Z bb)\^ + j d {PS,) 5^|Acai(u - W/Z hi + A;)|2 

= ^' + ^. (2.5) 

where the term integrated over the phase space measure d {PS3) corresponds to the virtual 
one-loop corrections (three final-state particles), while the one integrated over the phase space 
measure d {PS4) corresponds to the real tree level corrections with one additional emitted 
parton (four final-state particles). The sum ^ indicates that the corresponding amplitudes 
squared, | Airt{reai)(^j W/Z bb{+k))\'^, have been averaged over the initial-state degrees of 
freedom and summed over the final-state ones. The final phase space integration have been 
performed using Monte Carlo techniques using the adaptive multi-dimensional integration 
routine VEGAS [71] • 

Intermediate stages of the calculation of Safj^'~' contain both ultraviolet (UV) and infrared 
(IR) divergences, whose calculation will be discussed in detail in Sections 12.2.21 12. 2. 3^ 12.31 
(for Wbb production) and 12.41 (for Zbb production). 

2.2.1 Renormalization Scale Dependence of the NLO Cross Sec- 
tion 

We observe that the scale dependence of the total cross section at NLO is dictated by 
renormalization group arguments. In order to assure the renormalization scale independence 
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of the total cross section at 0{a^), f^^'~'{xi, X2,/i) has to be of the form: 

/|[LO(a:i,X2,/i) = fl,{x^,X2) + ^;(a:i,X2)ln (^^^ , (2.6) 

where, taking into account that other sources of renormahzation scale dependence in Eq. 02.41) 
are ^^(yu) and the PDFs Tf'^lx, fi), we can prove that: 

f!^ix,,x,) = 2\A7rbofl^''{xuX2)-Y,\ f dz,Pik{zi)fj:^{x,z,,x,) 



k '-^p 

1 

dZ2Pjk{z2)flf{Xl, X2Z2) 



(2.7) 



where p = (2mf, + MvY / s {V = W, Z), Pij are the Altarelli-Parisi splitting functions as 
presented in Eqs. fl2^ . ([Mil), ^12188]) and IK92^ . and bo is determined by the one- loop 
renormahzation group evolution of the strong coupling constant a^: 



-boal + 0{al), bo = j^(^^N-^nify (2.^ 



(iln(/i2) 

where = 3 is the number of colors and riif = 5 the number of light flavors. To write 
Eq. ( 12. 7p we have assumed that the 6-quark mass does not run, given the mild dependence 
of the cross section on it. The origin of the rest of terms in Eq. ( 12. 7p will become clear in 
Sections 12.31 and 12.41 where we describe in detail the calculation of virtual, real and total 
cross section corrections for both Wbb and Zbb production respectively. 

2.2.2 Calculating the Virtual Cross Section o"^^''* 

The 0{as) virtual corrections to the hadronic W/Z bb production tree level processes consist 
of self-energy, vertex, box and pentagon diagrams. The contributions to a^"'* in Eq. (12.51) 
can be written as: 

^1 Airt(^J ^ W/Z bb)\' = (A^L + AIAd) = E E^^^ {^oA) , (2-9) 

D D 

where Aq is the tree level amplitude and Ad denotes the amplitude for the one- loop diagram 
D, with D running over all self-energy, vertex, box and pentagon diagrams corresponding to 
the ij-initiated subprocess. 

The amplitude for each virtual diagram (Ad) is calculated as a linear combination of 
Dirac structures with coefficients that depend on both tensor and scalar one-loop Feynman 



12 



integrals with up to five denominators. We solve the one-loop integrals in the coefficients 
either at the level of the amplitude or at the level of the amplitude squared (see Eq. fl2.9p ). 
These two independent approaches allow us to thoroughly cross check the calculation of each 
individual diagram. Indeed, the tensor structures present in the one-loop integrals of the 
amplitude are typically different from the ones present in the amplitude squared, as one 
can perform non-trivial reductions of the latter by canceling dot products of the integration 
momentum in the numerator with denominators in the Feynman integrals. In this way, 
the final analytical expression of a given diagram ends up being represented in terms of 
different building blocks. A possible incorrect relation between the building blocks would 
then naturally produce a discrepancy between the two approaches. 

Tensor and scalar one-loop integrals are treated as follows. Using the Passarino-Veltman 
(PV) method [721 [73], the tensor integrals are expressed as a linear combination of tensor 
structures and coefficients, where the tensor structures depend on the external momenta and 
the metric tensor, while the coefficients depend on scalar integrals, kinematics invariants 
and the dimension of the integral (for a more detailed description of the technique see 
Appendix [D]). Numerical stability issues may arise at this level as a consequence of the 
proportionality of the tensor integral coefficients to powers of inverse Gram Determinants 
(GDs), specially when considering a full set of independent momenta {pa,} (i = 1, . . . , 4) of 
the ij — > W/Z bb phase space, which is defined by GD = det (p^. ■ Pa^)- The problem becomes 
more serious for higher rank tensor integrals, since the higher the rank of the original tensor 
integral, the higher the inverse power of the GD that appears in the coefficients of its tensor 
decomposition. 

To illustrate the problem, we parametrize the GD appearing in pentagon tensor integrals 
in terms of the W/Z bb phase space variables as 

GD = - ~ ^'^'^^^ ^^^'^ m + hi)' - ^M'y{s + -s,i)] s -s,i sin' 6,, sm' <f>,-, sin' 9 , 

(2.10) 

where s = xiX2Sh is the partonic center-of-mass energy squared. My is the mass of the weak 
vector boson {V = W,Z), and the W/Z bb phase space has been expressed in terms of a 
time-like invariant 5^5= {pb + Piy, polar angles {0,6^1) and azimuthal angles (0,0^5) in the 
center-of-mass frames of the incoming partons and of the bb pair, respectively. As can be seen 
in Eq. fl2.10p . the GD vanishes when the set of momenta become degenerate or co-planar, 
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for example at the boundaries of phase space. Near these regions of phase space it can 
become arbitrarily small, giving rise to spurious divergences which cause serious numerical 
difficulties, since large cancellations then appear between various parts of the calculation at 
the numerical level. The probability that the Monte Carlo integration hits a point close to 
these regions of phase space is not negligible and these points cannot just be discarded. 

The numerical instabilities we just discussed can be considered as "spurious" or "unphys- 
ical" divergences, since it is well known that only two-particle invariants can give rise to a 
physical singularity. Indeed, these spurious divergences cancel when large sets of diagrams 
are combined [M], such as, for example, when one combines gauge invariant sets of color 
amplitudes (i.e. amplitudes with a common color factor). Explicit analytic cancellations have 
been found for example when using helicity amplitudes and the helicity product formalism 
(see for example Refs. mainly because certain GD can be decomposed in terms of 

helicity products. As we have expressed our calculation in terms of kinematical invariants, 
the full cancellation only occurs between numerator and denominator at the numerical level, 
often between fairly large expressions. Nevertheless, when we consider gauge invariant sets 
of color amplitudes (as the ones presented in Tables 12.1112.51) and full analytical reductions of 
all tensor integrals, we find cancellation of some powers of GD, which improves the behavior 
of the numerical code so that we can integrate, using Monte Carlo techniques, over the entire 
phase space, to obtain statistical errors from the numerical integration below 0.1%. 

The fully reduced numerical amplitudes are often more demanding computationally, and 
because of that we have built numerical codes that use them only when close to regions of 
phase space where certain GD is small. With this the computer needs are reduced. All 
this was found particularly useful when considering E-PV functions (see Appendix |D]), and 
probably it would break down if one were to extend this technique to processes with even 
more legs, where probably using other techniques would be necessary. 

We also checked parts of our result by using unitarity techniques |5l], specifically the 
quadruple-cut technique [23]. As shown by Britto, Cachazo and Feng (BCF), from any set 
of Feynman diagrams (or more generally from any tensor integral [71]) one can extract the 
coefficient of a given scalar box integral by cutting the four corresponding propagators (see 
Fig. 12.11) . i.e. by replacing — m? + ie) 2'k5^^\P — ni^) for each cutted propagator of 
momentum i and mass m. This effectively freezes the momentum integration, and replaces 
it by a set of algebraic equations which determine the loop momentum entirely. We solved 
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this set of equations by using a BCF ansatz [23], and then compared the result to the 
corresponding box coefficient extracted from our analytic expression, and found agreement 
(for more details and specific solutions for the topology in Fig. 12.11 see Appendix [E]). This 
is a rather non-trivial check for the set of E-PV and D-PV functions (see Appendix [D]) 
we have employed at different stages, since they all contribute to the coefficients of the 
scalar D-functions occurring in the one- loop W/Z bb amplitudes. For instance, it has been 
particularly useful in the case of box diagrams like the one shown in Fig. 12.11 since this 
diagram and related ones contain up to DA-PV functions that cannot be reduced even at the 
level of the amplitude squared. Since they involve up to four powers of inverse GDs, they 
are particularly subject to numerical instabilities and it is important to have their analytic 
expressions as compact as possible. 




Figure 2.1: Quadruple cut [23] check of the calculation of a box diagram involving a top- 
quark loop. It corresponds to two Feynman diagrams {b[^} in Fig. I2.18P given by the two 
possible orientations of the fermion line. 



After the tensor integral reduction is performed, the fundamental building blocks are 
one-loop scalar integrals with up to five denominators. They may be finite or contain both 
ultraviolet (UV) and infrared (IR) divergences. 

The UV singularities of the virtual cross section are regularized in d = 4 — 2euv 
dimensions and renormalized by introducing a suitable set of counterterms, while the residual 
renormalization scale dependence is checked from first principles using renormalization group 
arguments as in Eq. 02.61) . The IR singularities of the virtual cross section are extracted in 
d = 4 — 2eiPi, dimensions and are canceled by analogous singularities in the C(«s) real cross 
section. 
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In our calculation we treat 75 according to the naive dimensional regularization approach, 
i.e. we enforce the fact that 75 anticommutes with all other 7 matrices in d = 4 — 2e 
dimensions. This is known to give rise to inconsistencies when, at the same time, the d- 
dimensional trace of four 7 matrices and one 75 is forced to be non-zero (as in d = 4, where 
Tr{'~f^'~f'^'~f'''~f'^'~f^) = Aie^'^^'^) [7S|. In our calculation, both UV and IR divergences are handled 
in such a way that we never have to enforce simultaneously these two properties of the Dirac 
algebra in d dimensions. For instance, the UV divergences are extracted and canceled at 
the amplitude level, after which the c? — > 4 limit is taken and the renormalized amplitude is 
squared using d = 4. Thus, all fermion traces appearing at this point are computed in four 
dimensions and therefore have no ambiguities. 

The finite scalar integrals are evaluated using the method described in Ref. [73] and cross 
checked with the numerical package FF [57]. The singular scalar integrals are calculated 
analytically by using dimensional regularization in d = 4: — 2e dimensions. The most difficult 
integrals arise from IR divergent pentagon diagrams with several external and internal 
massive particles. We calculate them as linear combinations of box integrals using the 
method of Ref. [76l [77] and of Ref. [73]. Details of the box scalar integrals, including two 
that we calculated explicitly since they were not previously in the literature, and of the 
pentagon reduction used in this calculation are given in Appendix [Cl as well as the set of 
UV- and IR-divergent three and two point functions. 

We note that the tree level amplitude ^0 in Eq. (12. 9p has generically to be considered 
as a (i-dimensional tree level amplitude. This matters when the Ad amplitudes in Eq. (12.91) 
are UV or IR divergent. Actually, as will be shown in the following, both UV and IR 
divergences are always proportional to pieces of the tree level amplitudes and they can be 
formally canceled without having to explicitly specify the dimensionality of the tree level 
amplitude itself. After UV and IR singularities have been canceled, everything is calculated 
in d = 4 dimensions. 

2.2.3 Calculating the Real Cross Section d-^^^^ 

The NLO real cross section a'^^^^ in Eq. (12. 5p corresponds to the 0{as) corrections to 
ij — i> W/Z bb due to the emission of an additional real extra parton, i.e. to the process 
ij — > W/Z bb + k. It contains IR singularities which cancel the analogous singularities 
present in the 0{as) virtual corrections and in the NLO PDFs. These singularities can be 
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either soft, when the emitted extra parton is a gluon and its energy becomes very small, or 
coUinear, when the final state parton is emitted collinear to one of particles in the initial 
state. There is no collinear singularity from the final b and b quarks, because the 6-quark 
mass regularize the collinear divergence. 

These IR singularities can be conveniently isolated by slicing the phase space of the final 
state particles into different regions defined by suitable cutoffs, a method which goes under 
the general name of Phase Space Slicing (PSS). The dependence on the arbitrary cutoffs 
introduced in slicing the final state phase space is not physical, and cancels at the level of 
the total real hadronic cross section, i.e. in a^^^^, as well as at the level of the real cross 
section for each separate subprocess. This cancellation constitutes an important check of 
the calculation and will be discussed in detail in Sections 12.3.31 and 12.4.51 

We have calculated the cross section for the processes 

i{qi) + j{q2) b{pb) + biPi) + V{py) + h{k) , 

with qi + q2 = Pb+Pb~^Pv + k, using the two-cutoff' PSS method, which includes cutoffs of the 
soft and collinear kind. This implementation of the PSS method was originally developed 
to study QCD corrections to dihadron production [TS] and has since then been applied to a 
variety of processes. (A nice review about it can be found in Ref . [71] , to which we refer for 
more extensive references and details.) 

In the following, we briefly review the structure of the real calculation using the two-cutoff 
PSS. We mention that the soft and collinear kernels employed throughout our calculation 
have been used also in the NLO calculation of ttH production at hadron colliders, where 
results have been checked using a PSS method with one-cutoff and a dipole cancellation 
method [801 El ESI [83]. Although the processes we are considering are different, the 
kinematics are equivalent, and the color structure and IR behavior are the same, so 
necessarily their soft and collinear kernels are the same. 

Phase Space Slicing method with two cutoffs 

The general implementation of the PSS method using two cutoffs proceeds in two steps. 
First, to isolate the soft singularities of a final state extra gluon we introduce an arbitrarily 
small soft cutoff 6s and we separate the overall integration over the phase space of that gluon 
into two regions according to whether the energy of the final state gluon {k^ = Eg) is soft, 
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i.e. Eg < Ss\^/2, or hard, i.e. Eg > Ssy/s/2. The partonic real cross section of Eq. fl2.5l) can 
then be written as: 

±.rcal ±.soft I '-hard In ^ ^\ 

(T = cr + (7 , \'^-^^) 

where obtained by integrating over the soft region of the gluon phase space, and 

contains all the IR soft divergences of a'^^^^. To isolate the remaining collinear divergences 
from a^'^'^'^, we further split the integration over the phase space of any final state parton 
according to whether the parton is [a^'^^d/coii-^ -g j^Q^ (^^hard/r,on-coii-^ emitted within an angle 
9 from the initial state partons such that (1 — cos^) < 5c-, for an arbitrary small collinear 
cutoff Sc- 

^hard ^hard/coll _j_ ^hard/non—coll ^2 \2i^ 

The hard non-coUinear part of the real cross section, (J ^ ^ IS finite and can be 

computed numerically, using standard Monte Carlo techniques. 

On the other hand, in the soft and collinear regions the integration over the phase space 
of the emitted gluon or quark can be performed analytically, thus allowing us to isolate and 
extract the IR divergences of 0"''°^'. More details on the calculation of a^°^^ and a^"-'^'^ for each 
relevant subprocess will be given in Sections 12.3.21 and 12.4.41 for the hadronic production of 
Whh and Zhh respectively. 

The cross sections describing soft, collinear and IR-finite radiation depend on the two 
arbitrary parameters 5^ and 5c- However, in the total real hadronic cross section (j^^^^, after 
mass factorization, the dependence on these arbitrary cutoffs vanishes for sufficiently small 
values of the cutoffs. 
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Figure 2.2: Tree level Feynman diagrams for gg' Whh. 
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2.3 Calculation of NLO QCD Corrections to Wbb 
Production at Hadron Colliders 



In this Section we present in detail the calculation of the partonic total or differential cross 
section a^^^{ij Wbb) [8l|, as defined in Eqs. ([23D and fl23D . 

There is only one subprocess contributing to Wbb at LO, qq' Wbb, where q and 
q' represent quarks or antiquarks of up-type and down-type, respectively. We neglect 
contributions from third generation initial quarks as they are suppressed by either the initial 
state quark densities (PDFs) or by the corresponding Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements. 

The contributing tree level Feynman diagrams are shown in Figure 12. 2[ Given the 
assignment of momenta 

q{qi)q'{q2) b{ph) + Kp-,) + w{p^-) , 

the LO amplitude can be written as 

Aoiqq' - Wbb) = ^gl^V,,. e;{p„) J^^^2 ^^^'^i 

{-q2 + Pw)^ 

+%7'^ /^~^" 7"(1-75K1 , (2-13) 

[qi - Pw)^ J 

where Qs and Qw are the strong and weak coupling constants, respectively, = are 
given in terms of the Gell-Mann matrices A" and Vqq' are the entries of the CKM mixing 
matrix. (For more details see Appendix iBl) 

The partonic LO cross section is obtained by integrating |^oP over the Wbb final state 
phase space: 

aLO = I d{PS,)Y,\M\ (2.14) 

where the sums indicates average over initial and sum over final spins and colors of the 
fermion lines, as well as sum over polarizations of the vector boson. As we are considering 
an on-shell gauge boson, we have summed over its polarizations according to the prescription 
used for massive vector bosons, i.e.: 

E^m(^)^^(^) = -^/^^ + ^' (2.15) 
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where is the mass of the V weak boson {V = W, Z). For Wbb production the second 
term in Eq. fl2.15l) does not contribute, because in our calculation the W boson only couples 
to the initial massless fermion line. The full expression has to be considered, however, when 
calculating Zbb production. 

At NLO one has to consider three processes: qq' Wbb, which contributes both at 
0{a'j.) and at 0{a^) through the one-loop 0{as) virtual corrections, and the real 0{as) 
corrections, due to qq' Wbb + g and q{q)g Wbb + q'{q') which contribute at 0{a^). In 
the following sections we will discuss in detail the structure of both virtual and real 0{as) 
corrections. 




Figure 2.3: Gluon (S*} ' ) and quark (5*2 ' ) self-energy corrections contributing to the 

qq' — >• Wbb subprocess at NLO. The shaded blobs denote standard one- loop QCD corrections 
to the gluon and quark propagators, respectively. 



2.3.1 Virtual Corrections to qq' Wbb 

The 0{as) virtual corrections to the qq' — > Wbb tree level process consist of the self-energy, 
vertex, box and pentagon one-loop diagrams illustrated in Figures 12.31 12. 4[ 12.51 and 12.61 The 
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Figure 2.4: 0{as) vertex corrections contributing to the subprocess qq' — > Wbb at NLO. The 
shaded blobs denote standard one- loop QCD corrections to the qqg, bbg and qq'W vertices, 
respectively. 



contributions to the virtual amplitude squared of Eq. (12. 5p can then be written as: 

^1 Ah.t(gg ^ wbb)\' = E E i^^-^D + A^d) = $^^27^e (a^I,) , (2.I6) 

D D 

where Aq is the tree level amplitude given in Eq. fl2.13p and corresponding to the diagrams 
shown in Figure fI72[ and Ad denotes the amplitude for the one-loop diagram D, with D 
running over all self-energy, vertex, box and pentagon diagrams illustrated in Figures 12. 3^ 
[2ai[23]and[2Sl 

The amplitude of each virtual diagram {Ad) is calculated as described in Section [2.2.21 
Inserting all diagram contributions into Eq. (12.161) . we obtain the complete 0{al) 
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Figure 2.5: 0{as) box diagram corrections contributing to the qq' — > Wbb process at NLO. 
The term Crossed refers to the box diagrams, -B^'^^ obtained from the -83^'^"* boxes by flipping 
the fe-quark fermion hne. 




Figure 2.6: 0{as) pentagon diagram corrections contributing to the qq' Wbb process at 
NLO. 



contribution to the virtual amplitude squared, and integrating over the final state phase 
space we calculate a'^^^ in Eq. (12. 5p . 

Results for the renormalization of the one-loop corrections are shown in Section 12.3.1.11 
The structure of the IR singular part of the virtual cross section is presented in Section [2.3.1.2l 
while the IR singularities of the real cross section are discussed in Section 12.3.21 The 
explicit cancellation of IR singularities in the total inclusive NLO cross section is outlined 
in Sections 12.3.21 and 12.3.31 
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2.3.1.1 Virtual corrections: UV singularities and counterterms 



The UV singularities of the 0{al) total cross section originate from the self-energy and 
vertex virtual corrections shown in Figures [23] and [231 These singularities are renormalized 
by introducing counterterms for the wave function of the external fields {6Z2''\ ^Z^"*) and 
the strong coupling constant ((5Zq,J. If we denote by Auv(5'j ' ) and Auv(V^ ' ) the UV- 



divergent contribution of each self-energy (S'l^'^'*) or vertex diagram iyr^') to the virtual 
amplitude squared (see Eq. fl2.16p ). we can write the UV-singular part of the total virtual 
amplitude squared as: 

2 3 



^{l,2)^ 



UV 1 2 
virt I 



Y)A,? Y. ^"v(S.'" + Sf ) + Auv(f <" + vP) (2.17) 
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We denote by \Ao\'^ the matrix element squared of the tree- level amplitude for qq' 
computed in d = 4 dimensions (see Eq. f l2.13p and Section r2.2.2p . 

The UV-divergent contributions due to the individual diagrams are given by: 
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(2.18) 



where iV = 3 is the number of colors, riif = 5 is the number of light flavors and AT^ and A/f, 
are standard normalization factors defined as: 



/ 4:71 
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V ^6 



Til + e) 



(2.19) 



Moreover, we define the required counterterms according to the following convention. For 
the external fields, we fix the wave-function renormalization constants of the external fields 
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(Z2^ = 1 + SZ2 \ i = q,b) using on-shell subtraction, i.e.: 



6Z. 



aA /N 1 \ / 1 



uv \4nJ \2 2N J \e 



-vv 



We notice that both 6z!f^ and 6Z2 \ as well as some of the vertex corrections {Vi^''^^ 

(1 2) 

and V2 ' ), have also IR singularities. In this section we limit the discussion to the UV 
singularities only, while the IR structure of this counterterm will be included in the IR- 
singularities shown in Section [2.3.1.21 



Finally, for the renormalization of we use the MS scheme, modified to decouple the 



top quark |85]. The first ri// light flavors are subtracted using the MS scheme, while the 
divergences associated with the top-quark loop are subtracted at zero momentum: 



:2.2r 



such that, in this scheme, the renormalized strong coupling constant as evolves with riif = 5 
light flavors, as justified by the energy scale of the processes under consideration. 

It is easy to verify that the sum of all the UV-singular contributions as given in Eq. (12.171) 
is finite. We also notice that the left over renormalization scale dependence, due to the 
mismatch between the renormalization scale dependence of Auv(>S'i) and 6{ZaJ, is given by: 

SAP^(-|....^).n(^) , (...) 

and corresponds exactly to the first term of Eq. (12.71) . as predicted by renormalization group 
arguments. We note that the presence of s in the argument of the logarithm of Eq. (12.221) 
has no particular relevance. Choosing a different argument would amount to reabsorbing 
some //-independent logarithms in /^"' of Eq. (12.61) . 

2.3.1.2 IR singularities 

This section describes the structure of the IR singularities originating from the 0{as) virtual 
corrections. The virtual IR singularities come from the following set of diagrams: vertex 
diagrams V^/^'^^ and V2^''^\ box diagrams B^''^\ B2, B^'"^^ and B^^''^^ and pentagon diagrams 
Pi and Pic, and from the wave function renormalization of the external fields, Sz!f^ and 
5Z2^ . After grouping all IR poles from this diagrams we obtain the total structure of the IR 
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singularity of the one- loop virtual corrections to qq' Wbb. Before writing such expressions, 
let us introduce the following set of kinematical variables: 

(gi + g2)' , 

- (91 -Pbf = 2qi-pb , 

'ml - {(12 -Pbf = 2q2-pb , 
^b - (Qi -Pbf = 2qi-p-b , 



s 

T2 
T3 

^bb 



{Pb+Pbf = 2pb-p-b + 2ml . 



Summing all the IR-divergent contributions yields: 



5:i-^*.i^=(^)^'5:i-^oi 



^virt 



with 



N 



(2.23) 



^2.24) 
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^bb 
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rat 
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^bbf^l 



A 



bb 
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riT2 



where we have used the kinematical invariants presented in Eq. (12.231) while Phi and are 
defined by: 



Pbb 
A 



bb 




(2.26) 



and (5^?;^ is a finite term that comes from having factored out a common factor A/b and is 
given by: 

Ct=fiV-il [^nf4')l+l[lln^f4')l . (2.27) 



A^ 



mt 



mt 



Before finishing let us write the IR-pole contributions to the counterterms SZ2''^ and 



SZ. 



(2.28) 
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ywvw V 



q,q>q',q' 



Figure 2.7: 0{as) real corrections: examples of initial and final real gluon emission and q{q)g 
initiated subprocess. 



In Sec. 12.3.21 we will show how the IR singularities of the real cross section exactly cancel 
the IR poles of the virtual cross section (see Eqs. f l2.35l) -( l2l36l) ). as predicted by the Bloch- 
Nordsieck [Sn] and Kinoshita-Lee-Nauenberg [HZl [SB] theorems. 

2.3.2 Real Corrections to Wbb Production 

The 0{as) corrections to qq' —* Wbb due to real gluon emission (see Figure 12. 7p give rise 
to IR singularities which cancel exactly the analogous singularities present in the 0{as) 
virtual corrections (see Sec. I2.3.1.2p . We also have real contributions from the subprocess 
q{(1)9 ~^ Wbb + q'{q') that give rise to IR singularities of the coUinear kind. We present 
results for the latter at the end of this section. 

We have calculated the cross section for the process 

q{qi) + q'iq2) ^ b{pt,) + b{pi) + W{p^) + g{k) (2.29) 

using the two-cutoff PSS method, as presented in Section [2.2.31 In the following subsections 
we explain in detail how we have applied this method to the calculation of the real 
contributions to Wbb hadronic production. We will present details of the calculation of 
the pieces of (7^°^"^ introduced in Section 12.2.31 namely a^°^^ and a^"-^'^. 
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Soft gluon emission 



The soft region of the qq' — > Whh + g phase space is defined by requiring that the energy of 
the gluon satisfies: 

Eg < , (2.30) 

for an arbitrary small value of the soft cutoff Sg- In the limit when the energy of the gluon 
becomes small, i.e. in the soft limit, the matrix element squared for the real gluon emission, 



real 



assumes a very simple form, i.e. it factorizes into the LO matrix element squared 



times an eikonal factor $ 



eik- 



where the eikonal factor is given by: 



eik 



N 
+ 



2 

mt 



+ 



+ 



T2 



{pi-kf {qi-k){pb-k) {q2-k){pi-k)_ 



(2.31) 



(2.32) 



2N 
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+ 



mr 



''>bb 



2ml 



{pb-ky {Pb-ky {qrk){q2-k) {pb-k){p-b-k) 



T2 



{(li-k){pb-k) {qi-k){p-b-k) {q2-k){pb-k) {q2-k){p-b-k) 
where we have used the kinematical invariants defined in Eq. (12.231) . Moreover, in the soft 
region the qq' —>■ Whh + g phase space also factorizes as: 



d{PSi){qq' Whh + g) 



soft 



d{PSs){qq' ^Whb)d{PSg)soft 
d{PS,){qq' - W^fe&)77T37;iZTv77^^(5.^ " Eg) 



(2.33) 



(27r)('^-i)2Eg 

where d{PSg)soft denotes the integration over the phase space of the soft gluon. The parton 
level soft cross section can then be written as: 

a-/* = (47ra,)/i2^ j d{PSs)J2\Ao\' J d{PSg)soft^,ik ■ (2.34) 

Since the contribution of the soft gluon is now completely factorized, we can perform the 
integration over d{PSg)soft in Eq- (12.341) analytically, and extract the soft poles that will 
cancel X™^ and X™* of Eq. (I2.25p . The integration over the gluon phase space in Eq. (12.341) 
can be performed using standard techniques and we refer to Refs. [791 [89] for more details. 
For sake of completeness, in Appendix [F] we give explicit results for the soft integrals used 
in our calculation. 
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Finally, the soft gluon contribution to o'gl'^^ can be written as follows: 



soft 



where e corresponds to eiR of Eq. fl2.24p and 
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(2.36) 



while A/fe is defined in Eq. (12.191) . Li2 denotes the dilogarithm as described in Ref. [90] and 
Xri* and X™^ are given in Eq. (I2.25p . We have used the kinematical invariants defined in 
Eq. (I2:23|) . (3bb and A^s are defined in Eq. ^2^, 

%b 



abb 



l+/5b-b "^bKb-1) 



1-A 



and 



bb 



'bb 



2(ab-bPi;-p-? 



(2.37) 
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while, for any initial parton i and final parton /, the function Fij can be written as: 

F., = Viln^fi±|t) (2.38) 

+2L J-Ml^i^^V 2Li/ '3/(1 +»^»-/) 



1-/9/ J \ l — PfCos6if 

where 6if is the angle between partons i and / in the center-of-mass frame of the initial state 
partons, and 

I ~ 

All the quantities in Eq. (12.380 can be expressed in terms of kinematic invariants, for details 
see Appendix [Fl 

As can be easily seen from Eqs. fl2.25p and fl2.36p . the IR poles of the virtual corrections 
are exactly canceled by the corresponding singularities in the soft gluon contribution. The 
remaining IR poles in a^°f^ will be canceled by the PDF counterterms as described in detail 
in Sec. 12^:31 

Hard gluon emission 

The hard region of the gluon phase space is defined by requiring that the energy of the 
emitted gluon is above a given threshold. As we discussed earlier this is expressed by the 
condition that 

E,>Ss^ , (2.40) 

for an arbitrary small soft cutoff 6s, which automatically assures that a'^"-^'^ does not contain 
soft singularities. However, a hard gluon can still yield singularities when it is emitted 
at a small angle, i.e. collinear, to a massless incoming or outgoing parton. In order to 
isolate these divergences and compute them analytically, we divide the hard region of the 
qq' — s> Wbb+g phase space into hard/ collinear and hard/non- collinear regions, by introducing 
a small collinear cutoff 5c- The hard/non- collinear region is defined by the conditions 

^>5. and ^>4. (2.41) 

The contribution from the hard/non- collinear region, (j-'^o.rd/non-coii ^ finite and we compute 
it numerically using standard Monte Carlo integration techniques. 
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In the hard/collinear region, one of the conditions in Eq. fl2.41l) is not satisfied and the 
hard gluon is emitted colhnear to one of the incoming partons. In this region, the initial-state 
parton i {i = q, q) is considered to split into a hard parton i' and a collinear gluon g, i ^ i'g, 
with = zpi and k = {l — z)pi. The matrix element squared for ij — * Wbb + g factorizes into 
the LO matrix element squared and the Altarelli-Parisi splitting function for i — > i'g, i.e.: 

.2R,Az) 



X]l-^real(U ^ Wbb + g] 



12 collinear 



{47ra.)Y,Y.\M^J^Wbb)\' 



z s 



(2.42) 



«9 



with Sig = 2pi-k, and Pn'^z) is the unregulated Altarelli-Parisi splitting function ioi q + g 
at lowest order, including terms of 0{e) as given by Pqq{z) = P^^ + eP^^ with 

1 + ^2 



Pqq{z) 



l-z' 
-Cf{1-z) 



^2.43) 



where Cp = 1/2{N — 1/N). Moreover, in the collinear limit, the qq' —>■ Wbb + g phase space 
also factorizes as: 

d{PS,){zj Wbb + g) '"'^^^ d{PS,W3 - Wbb) e (Eg- 6^^^) x (2.44) 



(27r)(^-i)2E3 



^(cos^, 



«9 



r(l-e) (4.)'^_,^_^^_^„j_^)^^j-.,/(i^^,|_^_^ 



r(l - 2e) 167r2 

where the integration range for Sig in the collinear region is given in terms of the collinear 
cutoff, and we have defined s' = 2pii ■ pj. The integral over the collinear gluon degrees 
of freedom can then be performed separately, and this allows us to extract explicitly the 
collinear singularities of a^'""'^. (jhard/coii ^^^j^j^g Qf ^^j-^g form [THl 1^ : 

,2\ n 



hard/ coll 



r(l-e) /47r/i2 



_27rr(l - 2e) V 



dz 



[I - zf s' 
2z J 



(2.45) 



P,,{z)a'^''{i'j ^Wbb) + {i^ j) \ . 



The upper limit on the z integration ensures the exclusion of the soft gluon region. As usual, 
these initial-state collinear divergences are absorbed into the parton distribution functions 
as will be described in detail in the Sec. I2.3.3[ 
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The tree level processes (q, q)g — > Wbb + (q', q') 

The extraction of the collinear singularities of cf'^^^^ is done in the same way as described 
in the previous subsection for the qq' initial state. In the collinear region, cos 9iq > 1 — 5c-, 
the initial state parton i with momentum qi is considered to split into a hard parton i' and 
a collinear quark q, i ^ i'q, with gj/ = zqi and k = {1 — z)qi. The matrix element squared 
for ij Wbb + q factorizes into the unregulated Altarelli-Parisi splitting functions in d 
dimensions: Pu' = P^^, + eP/j, and the corresponding LO matrix elements squared. The 
ij — > Wbb + q phase space factorizes into the i'j Wbb phase space and the phase space 
of the collinear quark. As a result, after integrating over the phase space of the collinear 
quark, the collinear singularity of a^^g^^ can be extracted as: 



a. 



^^coll 
19 



27rr 1 -e 



-\ 1 5; 



dz 



z? s' 



2z 



mt 



X 



[P,,{z)a]^^^{q{q,)q{q,,)^Wbb)\ 



(2.46) 



The collinear radiation of an antiquark in qg Wbb + g is treated analogously. In the case 
of (g, q)g Wbb + (g, g) we have the possible splitting g qq. The 0(1) and 0(e) parts of 
the corresponding splitting function are: 



-z(l-z) 



(2.47) 



Again, these initial state collinear divergences are absorbed into the parton distribution 
functions as will be described in detail in Section 12.3.31 



2.3.3 Total Cross Section for pp{pp) Wbb 



As described in Sec. 12. 1^ the observable total cross section at NLO is obtained by con- 
voluting the NLO parton level cross section with the NLO parton distribution functions 
!F^'^{x, /i), thereby absorbing the remaining initial-state singularities of Sa^^^ into the quark 
distribution functions. This can be understood as follows. First the parton cross section is 
convoluted with the bare quark/antiquark distribution functions jF^'|(a;) and subsequently 
jF^'l(x) is replaced by the renormalized quark/antiquark distribution functions jF^'|(x, /i) 



defined in some subtraction scheme. Using the MS scheme, the scale-dependent NLO quark 
distribution functions are given in terms of jF^'^x) and the QCD NLO parton distribution 
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function counterterms [TH] as follows: 



(a) For the case where an initial state gluon g splits into a qq pair {g qq): 



27r 



r 1 



1 rf^ 

2; 



(2.48) 



where P^^ is defined in Eq. f l2.47p . This is relevant to process q{q)g Wbb + q'{q') 
when the gluon becomes coUinear with the final massless parton. 



(b) For the case oi q ^ qg splitting: 




(2.49) 



where the 0{as) terms in the previous equation are calculated from the 0{as) 
corrections to the q ^ qg splitting, in the two-cutoff PSS formalism, and Pqq{z) 
is the Altarelli-Parisi splitting function of Eq. fl2.43p . This is relevant to process 
qq' — > Whh + g when the final gluon goes soft or when one of the initial partons 
become collinear with the gluon. 

When convoluting the parton cross section with the renormalized quark/antiquark 
distribution functions of Eq. fl2.49p . the IR singular counterterm, that is the first term of 
the RHS of Eq. fl2.49|] . exactly cancels the remaining IR poles of a^^q"" + d'^°^* and 0-'*"''''^/'^°". 
Finally, the complete 0{a1) inclusive total cross section for pp{pp) 
factorization scheme can be written as follows: 



Wbb in the MS 



a 



NLO I ^NLO 
Oqq' -+- Oqg+qg 



(2.50) 
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with a^^,° and a^g° defined in the following, Eqs. fl2.5ip and fl2.53|] respectively. 

: Yl [ dx,dx,J^^{Xu /i) J-J^^^(X2, fi) X2, fi) + <3-;f (Xi, X2, /i) + a' X2, /i)] 



NLO 
^qq' 
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(2.52) 



1 ^ 2) 
Xi,X2,/i) + (1 ^ 2)} . 



(2.53) 



We note that 0"^^° is finite, since, after mass factorization, both soft and coUinear 

and a^'^-^d/cou^ ^^^^ ^YiaX the second 



singularities have been canceled between + a' 
term in Eq. fl2.5ip . which is proportional to In {^'^^ corresponds exactly to the second and 
third terms of Eq. (12.71) . as predicted by renormalization group arguments. 

To finish this Section, and before we discuss in detail in Chapter [3] the numerical results 
for the NLO total cross section for pp — ^ Whh, we first demonstrate that cr^^'-' does not depend 
on the arbitrary cutoffs of the PSS method, i.e. on the soft and hard/collinear cutoffs 5s and 
5c- We note that the cancellation of the cutoff dependence at the level of the total NLO cross 
section is a very delicate issue, since it involves both analytical and numerical contributions. 
It is crucial to study the behavior of 0""^^° in a region where the cutoffs are small enough 
to justify the approximations used in the analytical calculation of the IR-divergent part of 
(jgl'^', but not so small to cause large numerical cancellations. The Monte Carlo phase space 
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Figure 2.8: Dependence of (T^^°(pp — * WW) on the 5^ PSS parameter, when be is fixed at 
be = 10"^. In the upper window we illustrate separately the cutoff dependence of the soft 
and hard-coUinear part (2 — 3, red dashed curve) and of the hard non-coUinear part (2 — >■ 4, 
blue dotted curve) of the real corrections to the total cross section. The 2 — > 3 curve also 
includes those parts of the 2 — >• 3 NLO cross section that do not depend on 6c and 6s, i.e. the 
tree level and one-loop virtual contributions. The sum of all the contributions corresponds 
to the black solid line. The lower window shows a blow-up of the black solid line in the 
upper plot, to illustrate the stability of the result. The error bars indicate the statistical 
uncertainty of the Monte Carlo integration. 



integration has been performed using the adaptive multi-dimensional integration routine 
VEGAS fTT]. 

Figures [2^8] and [2^9] illustrate the dependence of the total cross section on the two-cutoffs 
of the PSS method, using the setup outlined in Section 13.11 In Figure 12. 8[ we show the 
dependence of cr^^'-' on the soft cutoff, 6s, for a fixed value of the hard/coUinear cutoff, 
5c = 10"^. In Figure [2l9| we show the dependence of a^^° on the hard/collinear cutoff, 6c, for 
a fixed value of the soft cutoff, 6s = 10~^. In the upper window of Figure [278112. 9p we illustrate 
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Figure 2.9: Dependence of a^^°{pp — * Wbb) on the 6c PSS parameter, when 6s is fixed at 
6s = 10^^. In the upper window we illustrate separately the cutoff dependence of the soft 
and hard-coUinear part (2 3, red dashed curve) and of the hard non-collinear part (2 — > 4, 
blue dotted curve) of the real corrections to the total cross section. The 2 — > 3 curve also 
includes those parts of the 2 — >• 3 NLO cross section that do not depend on 6c and 6s, i.e. the 
tree level and one-loop virtual contributions. The sum of all the contributions corresponds 
to the black solid line. The lower window shows a blow-up of the black solid line in the 
upper plot, to illustrate the stability of the result. The error bars indicate the statistical 
uncertainty of the Monte Carlo integration. 



the cancellation of the 5^(4) dependence between asoft + (^hard/cou and ahard/non-cou, while 
in the lower window we show, on a larger scale, a^^° with the statistical errors from the 
Monte Carlo integration. As before, a^^° also includes the contribution from the LO and 
the virtual cross sections, which are both cutoff-independent. For 6s in the range 10"^ — 10^^ 
and 6c in the range 10~^ — 10~^, a clear plateau is reached and the NLO total cross section is 
independent of the arbitrary cutoffs of the two-cutoff PSS method. All the results presented 
in Chapter [3] are obtained using the two-cutoff PSS method with 6s = 10~^ and 6c = 10"^. 
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2.4 Calculation of NLO QCD Corrections to Zbb 
Production at Hadron Colliders 

In this Section we present in detail the calculation of the partonic total and differential cross 
section a^^'^lij —> Zbb) [22], which can be decomposed as in Eq. fl2.4l) . We work throughout 
in the 4-flavor number scheme, where only 4 massless quark flavors can be excited in the 
initial state, as we consider full 6-quark mass contributions to the partonic cross section. 

There are two subprocesses contributing to a^*^ in Eq. fl2.4p . namely qq Zbb and 
gg Zbb. 

The tree level Feynman diagrams contributing to the LO cross section for subprocess 
qq Zbb are shown in Figures 12.21 (with V = Z) and 12.101 when the Z weak boson is 
emitted from initial and final fermion lines respectively. In Figure [2. Ill we show the tree level 
Feynman diagrams contributing to the gg Zbb subprocess. The corresponding amplitudes 
are presented in Section 12.4.11 and more details are summarized in Appendix [Bl 




Figure 2.11: Tree level Feynman diagrams for gg Zbb. The circled crosses correspond to 
all possible insertions of the Z boson, each one representing a different diagram. 



We notice that the tree level diagrams for qq' — > Wbb are a subset of the diagrams for 
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qq Zbb, namely the ones in which the Z boson is emitted from the initial quarks legs. 
The same holds at NLO: the NLO QCD corrections to qq' Wbb discussed in Section fT^ 
are a subset of the NLO QCD corrections to qq Zbb discussed in this Section. The 
results derived for qq' Wbb can be automatically translated to the qq Zbb case by 
keeping the contributions of the vector and axial vector parts of the Wqq' vertex separate 
and substituting the corresponding values of vector and axial vector couplings of the Zqq 
vertex, as well as making the necessary changes to the PDFs of the initial partons in the 
hadronic cross section. 

At NLO one has to consider five subprocesses: qq —>■ Zbb and gg Zbb contributing at 
tree level to a^'~' in Eq. fl2.4p and at one-loop to a"''* in Eq. (12.51) . as well as qq — Zbb + g, 
Q{(l)g ~^ Zbb + q{q) and gg Zbb + g contributing at tree level to a''^**' in Eq. fl2.5p . 

The rest of this Section is organized as follows. We present in Section 12.4.11 results for 
the LO amplitudes. In Sections 12.4.21 and 12.4.31 we shall present results for the 0{as) virtual 
corrections to Zbb hadronic production. For details on the 0{a1) subprocess qq Zbb+g we 
refer the reader to Sections 12.3.21 and 12.3.31 as results are analogous to the 0{al) subprocess 
qq' Wbb + g. In Sections 12.4.41 and Section [2.4.51 we shall discuss the 0{al) subprocess 
gg Zbb + g. 0{al) real subprocesses initiated by q{q)g for Zbb are slightly different in 
structure with respect to the similar ones encountered in the Wbb, as the LO structure of 
the former is more complex than the latter. For that reason we discuss their structure in 
detail in Sections 12.4.41 and I2.4.5[ 

2.4.1 Tree Level Cross Section for Zbb Hadronic Production 

The contributing tree level Feynman diagrams for the LO qq Zbb process are shown in 
Figures 12.21 (with V = Z) and 12.101 for subprocess qq Zbb with the Z weak boson emitted 
from initial and final fermion lines, respectively. Given the assignment of momenta: 
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the LO amplitude can be written as: 
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(2.54) 



[{.-Pb-PzY - ml] 

where Qs and are the strong and weak couphng constant, respectively, = are 
given in terms of the Gell-Mann matrices A" and 9^ is the weak angle. The vector, gy, and 
axial, (7^4, couplings for the Zff vertex are given explicitly in Eq. (1A.28P (for more details 
see Appendices \M and |B]) . 



The tree level amplitude for the process 

9\qi) + 9\q2) ^ h{p,) + h{n) + Z{pz) , 

where gi + g2 = Pt+Pi+Pz and a, b denote the color of the incoming gluons, is obtained from 
the three classes of Feynman diagrams represented in Figure [2.111 identified as s— channel, 
t— channel, and u— channel diagrams, respectively. We find it convenient to organize the 
color structure of both the tree level amplitude and the one-loop virtual amplitude in terms 
of only two color factors, one symmetric and one antisymmetric in the color indices of the 
initial gluons. Following this prescription, the tree level amplitude for gg Zbb can be 
written as: 

A = ^r'[t",t1+v4;5H^",^'} , (2.55) 

where t"'''' = A"'^/2 are given in terms of the Gell-Mann matrices A'^''' §. A^^ and Aq""^ 

correspond to the terms in the amplitude that are proportional, respectively, to the abelian 

^We note that the one-loop virtual amplitude can be expressed in terms of the same antisymmetric color 
factor [i'',t''] and a symmetric color factor made of {t°-,t''} and J"''. 
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(or symmetric) and non-abelian (or antisymmetric) color factors and are given by: 

X = + ^o,u) , = A,s + l{Ao,t - Ao,u) , (2.56) 

where Aq^s, Ao^t, and Ao,u are the amphtudes corresponding to the sum of the s— channel, 
t— channel, and m— channel tree level diagrams in Figure [2. 11[ Explicit expressions for them 
are shown in Eqs. flB.3P - flB.5l) . 

Because of the orthogonality between symmetric and antisymmetric color factors, the 
tree level amplitude squared takes the very simple form: 



^{N^ - 1) (l^rf + lA^'l^) - ^{N^ - 1)|X|' 



(2.57) 



The partonic LO cross section is obtained by integrating |.4oP over the Zbb final state phase 
space, as shown in Eq. fl2.14p . We perform the sum over polarizations of the Z gauge boson 
using the prescription shown in Eq. fl2.15p . keeping both terms. 

When averaging over the polarization states of the initial gluons (in the gg Zbb case) , 
the polarization sum of the gluon polarization vectors, e^(gi,Ai) and €^(^2,^2), has to be 
performed in such a way that only the physical (transverse) polarization states of the gluons 
contribute to the matrix element squared. We adopt the general prescription: 

where i = l,2 and the arbitrary vectors rij have to satisfy the relations: 

< Yl e;.(g.,A,)e:(g^,A,) = , < ^ e^(g„ A,)e:(g„ A,) = , (2.59) 

Ai=l,2 Ai=l,2 

together with Ui ■ qj^O and rii^n-^. We choose ni = q-2 and n2 = qi, such that: 

^1^92;. + q2i,qiu 

Ai=l,2 

Finally, the entire calculation is performed using the Feynman gauge for both internal and 
external gluons (for more details see Appendix [A|) . 



Y: e,{q., A.)<(g., A.) = + ^^^^^^ ^ . (2.60) 

^ gi • 52 



2.4.2 Virtual Corrections to qq — > Zbb 

The 0{as) virtual corrections to the qq Zbb tree level subprocess consist of the self energy, 
vertex, box and pentagon diagrams illustrated in Figures when the Z weak boson is 
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Figure 2.12: Gluon (Sj'fO and 6-quark (S^'f ) ^'(a,) self energy corrections contributing 
to the qq Zbb subprocess at NLO, when the Z boson is emitted from the final fermion 
line. The shaded blobs denote standard one-loop QCD corrections to the gluon and quark 
propagators respectively. 



emitted from the initial fermion line (g or q) , and in Figures 12. 12112. 15| when it is emitted 
from the final fermion line {b or b). The contributions to the virtual amplitude squared of 
Eq. (12.51) can then be written as: 

^1 A.t(gg - Zbb)\' = J2T. (^oA + ^o^a) = E E2^^ (^oA) , (2.61) 

Di Di 

where Ao is the tree level amplitude, corresponding to the diagrams shown in Figures 12.21 
and 12.10] and Ad, denotes the amplitude for a one-loop diagram, with Di running over all 
self-energy, vertex, box and pentagon diagrams illustrated in Figures 12.31 12.41 12.51 12.61 12.12[ 
12031 EH and Em 

The calculation of each virtual diagram (Ad,) is performed in the way explained in 
Sections 12.2.21 and 12.3. 1[ As mentioned there, one of the most challenging parts of the 
calculation is related to controlling the spurious divergences that appear when reducing 
tensor Feynman integrals. This is specially true when considering pentagon diagrams, like the 
ones in Figures and [2.151 There were only two pentagon diagrams in the Wbb calculation, 
but that number is doubled just by considering the one-loop diagrams for qq Zbb, and 
we will see in Section 12.4.31 that when considering one-loop diagrams for gg Zbb one has 
to add twelve more pentagon diagrams, some of them containing up to EA-PV functions 
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Figure 2.13: 0{as) vertex corrections contributing to the qq Zbb subprocess at NLO when 
the Z boson is emitted from the final fermion fine or from a closed fermion line. The shaded 
blobs denote standard one-loop QCD corrections to the qqg (V^^V ), bbg {¥2.^' ) and qqZ 

(1 2) (1 2) 

(Vg.f' ) vertices respectively. V^.f' are b- and t- fermion loop vertices which are UV and IR 
finite (contributions of quarks from first and second family vanish). 



(see Appendix |D]). All this increases considerably the stability problem. In the case of 
pentagon diagrams it is convenient to reduce consistently all E-PV functions by canceling 
systematically, at the level of the amplitude squared in Eq. (12.611) . all possible vector products 
containing the loop momentum in the numerator with some denominators. This is possible 
as, in the pentagon topology of our process, each leg has an outgoing momentum which 
is on-shell, corresponding basically to one of the external initial or final particles of the 



41 




t!(l) r(2) 

Figure 2.14: 0{as) box diagram corrections contributing to the qq Zbb subprocess at 
NLO, when the Z boson emitted from the final fermion fines (6 or b). 



subprocess. One tfien ends witfi expressions for eacfi pentagon diagram containing purely 
scalar pentagon integrals, or tensor integrals witfi fewer tfian five denominators, improving 
considerably tfie befiavior of tfie numerical code. We cfiecked analytically tfiese reductions 
witfi tfie non-reduced expressions by using tfie full reduction of all tensor integrals to scalar 
integrals, and found agreement. Anotfier gain we fiad by doing tfie aforementioned reductions 
was a considerable speed up of all analytical and numerical computations. 

In Section [2.4.2.11 we present tfie UV singularity structure of a^"^{qq Zbb) as it fias a 
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Figure 2.15: 0{as) pentagon diagram corrections contributing to the qq 
at NLO, when the Z boson emitted from the final fermion lines (6 or b). 



Zbb subprocess 



slightly different structure to what is discussed in Sections 12 . 3 . 1 . 1 1 due to the Z boson emission 
from the final fermion lines. The IR divergent structure on the other side is analogous to 
that presented in Section [2.3.1.21 to which we refer the reader for details. 

2.4.2.1 Virtual corrections: UV singularities and counterterms 

The UV singularities of the 0{al) total cross section for qq — s> Zbb originate from self-energy 
and vertex virtual corrections shown in Figures 12.31 12.41 12.121 and 12.131 These singularities 
are renormalized by introducing counterterms for the wave function of the external fields 
{8Z^^ ^ 6Z^^\ the bottom-quark mass (^mb), and the strong coupling constant ((5ZqJ. If 
we denote by Auv(5'j ' ) and Auv(V^j ) the UV-divergent contribution of each self-energy 
{S\ ' or vertex diagram {V^ ) to the virtual amplitude squared (see Eq. (I2.6ip ). we can 
write the UV-singular part of the total virtual amplitude squared as: 

^K",y.r = El Ai^ (e a„v (sf ' + s.'^' + 4' + sg') 

+ X:Auv(v« + ^'nV.J' + v|') 

+ 2[(*^'')„v+(*^'')„v + *^-]} • 

We denote by l^loP the matrix element squared of the tree- level amplitude for qq — > Zbb, 
computed in d = 4 dimensions, as presented in Eq. (12.54p . 
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The UV-divergent contributions due to the individual diagrams are exphcitly given by: 
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(2.63) 



where Ms and Mh are standard normahzation factors defined in Eq. fl2.19p . 

We define the required counterterms according to the following convention. For the 
external fields, we fix the wave-function renormalization constants of the external fields 



1 + 5Z. 



2 5 



■q,b) using on-shell subtraction, given expressions as in Eq. fl2.20p . We 



notice that both ^Zg^'' and SZ2' , as well as some of the vertex corrections [Vi '''^' , V-!^.'',^' and 



(1.2) T.{1,2) 



/ 



V, 



(1,2) ..{l,2)^ 



have also IR singularities. In this section, we limit the discussion to the UV 



singularities only, while the IR structure of these terms are analogous to what is shown in 
Section EXH 

We define the subtraction condition for the bottom-quark mass nif, in such a way that 
is the pole mass, in which case the bottom-mass counterterm is given by: 
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rrib 
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(2.64) 



q(2) 
'^2;f 



this counterterm have been already included into the term Auv 
in Eq. (12:631) . 

The structure of the remaining counterterms that appear in this case is the same as 
presented in Section [2.3.1.11 and we refer the reader to it. 



2.4.3 Virtual Corrections to gg — > Zbb 

The 0{as) virtual corrections to the gg Zbb tree level subprocess consist of the self-energy, 
vertex, box, and pentagon diagrams illustrated in Figures [2. 16112. 191 The 0{a^) contribution 
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Figure 2.16: Gluon (5*1^/^) and quark (Sj^/^ S^,'s \ 5_|^(ff)^^) 0{as) self energy 

corrections contributing to the gg Zbb subprocess at NLO. The circled crosses correspond 
to all possible insertions of the Z boson, each one representing a different diagram. 



to the virtual amplitude squared of Eq. fl2.5p can then be written as: 
Y.lA.irtigg - Zbb)\' = (A^K,, + A^A.) = $^^27^e (A^K.) ' (2-65) 

where Aq is the tree level amplitude given in Eq. fl2.55p . while Ad^^ denotes the amplitude 
for a class of virtual diagrams that only differ by the insertion of the final state Z boson 
leg, i.e. Di^j = XIa: -^ij with Di = Si,Vi, Bi, Pi, j = s,t,u, and k running over all possible Z 
boson insertions, as illustrated in Figures [2.16112.191 

The calculation of each virtual diagram {Aoi^) is performed in the way we have explained 
in Sections I2.2.2[ 12.3.11 and 12.4.21 We refer to those sections for details. We note that the 
greatest complexity for virtual diagrams in our calculation is encountered in the gg Zbb 
subprocess. This is due to the higher number of one-loop diagrams that need to be considered 
and to the presence of two external gluons. For this reason we have chosen to organize the 
diagrams, at certain stages, into gauge invariant color amplitudes, that is, into coefficients 
of the same color structure. This allows a better handling of the spurious singularities 
and a natural way to make internal cross checks and cross checks with new techniques (see 
Section 12.2.21 and Appendix [E]) . 
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Figure 2.17: 0(0^) vertex corrections contributing to the gg Zbb subprocess at NLO. 
The shaded blobs denote standard one- loop QCD corrections to the ggg (y^J ), bbg J , 

'^r^tu)^ and V^^^^^^^) and bbZ {V^]^''^^ and V^^^^^^^) vertices. The circled crosses correspond to 
all possible insertions of the Z boson, each one representing a different diagram. 
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Figure 2.18: 0{as) box diagram corrections contributing to the gg Zbb subprocess at 
NLO. The circled crosses correspond to all possible insertions of the Z boson, each one 
representing a different diagram. 



We introduce the leading and sub-leading color factors: 

Cs = (l + ^) {N' - 1) , (2.66) 
and use them to group the different diagrams when interfered, as in Eq. fl2.65p . with the LO 
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Figure 2.19: pentagon diagram corrections contributing to the gg Zbb subprocess 

at NLO. The circled crosses correspond to all possible insertions of the Z boson, each one 
representing a different diagram. 



amplitudes in Eq. fl2.55p . Aq^ and Aq"''', independently. This leads to the seven sets of color 
amplitudes A'^''^"''^'' [i = 1, ... ,5), whose diagram content is shown in Tables 127X112.51 

Table 2.1: Diagram content of color amplitudes Af' and A""* whose color factor, when 
interfered with Aq^ or Aq"''^ respectively, is Ci. 



Color Factor: Ci 


self energy 


vertex 


box 


pentagon 




q{l,2) c.(l,2,3,4) 


..(1,2) ..(1,2,3) 

T/(l-2,3) 
9,{t,u) 


R r(1.2) 
^4,s,^(7,8),(t,„) 


Pl,it,u) 




C-(l,2) q(1,2) 

C.(l,2) c;,{l,2,3,4) 
'^3,{t,u)^ '-'■i,{t,u) 


..(1,2) ..(1,2) 

^1,S 5 ^2,S 5 

..(1,2) ..(1,2) 
•■^3,8 ' '■^6,5 ' 
..(1,2,3) ..(1,2,3) 
•^(7,8), (i.«)' ^9,it,u) 


-^3,55 -84,55 

o(l,2) 
^(7,8),(t,«) 


Pl,(t,u) 
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Table 2.2: Diagram of content color amplitudes A2 and A2 whose color factor, when 
interfered with or Aq""^ respectively, is C2. 



Color Factor: C2 


self energy 


vertex 


box 


pentagon 


Af 


C.(l,2) c;,(l,2,3,4) 


..(1,2,3) T.(l,2,3) 
''(7,8), (*,«)' ^9,it,u) 


p(l,2) p(l,2) 
-°(5,6),(t,«)' -°7,(t,-u)' 

5(1,2) 

-"(9,10), (t,«) 


-P(2,3,4),(t,«), 




0(1,2) 
<-'2 s 

^{1,2) 'c.(l,2,3,4) 
'^3,(t,-u)' '^4.,(t,u) 


T/(l,2) 
(2,3),s 
..(1,2,3) ..(1,2,3) 
'^(7,8),(t,-u)' ^9,{t,u) 


r(1,2) r(1,2) 
-"2,s 5 -0(5,6), (*,«)' 
r(1,2) 
9,(t,«) 


-P(2,3,4),(t,«) 



Table 2.3: Diagram content of color amplitude Af' whose color factor, when interfered with 
Aq'^ or ^q"-^ respectively, is C3. 



Color Factor: C3 


self energy 


vertex box pentagon 


At 


C.{1,2) c;,{l,2,3,4) 
'^3,(t,«)' '-'■i,{t,u) 


..(1,2,3) ..(1,2,3) ^(1,2) ^(1,2) p 
^(7,8),(t,u)^ ^9,it,u) -0(5,6),(t,«)' -0{9),(t,«) -f^(2,3,4),(t,«) 



Table 2.4: Diagram content of color amplitude A2°'^ whose color factor, when interfered with 
At or Aq''^ respectively, is NC2. 



Color Factor: NCo 



self energy 



vertex 



box pentagon 



j^nab 



(1,2) 



(1,2) 



(1,2) 



5,s 



T/(l,2) 
^4,(t,«)' 
T/(l,2,3,4) 
10,(t,«) 
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(1,2,3) 
t 
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Table 2.5: Diagram content of color amplitude whose color factor, when interfered with 
A^'' or A'f^"^ respectively, is (A^ - 2/N)C2. 



Color Factor: {N-2/N)C2 


self energy vertex 


box pentagon 


Af 


T/(l,2,3,4) 
W,{t,u) 


^(1,2,3) _ 



The structure of the UV singularities for the virtual cross section of the gg Zbb 
subprocess is presented in Section [2.4.3.1[ The structure of the IR singular part on the other 
hand is presented in Section 12.4.3.21 while the IR singularities of the real cross section are 
discussed in Section 12.4.41 The explicit cancellation of IR singularities in the total inclusive 
NLO cross section for gg Zbb is outlined in Sections 12.4.41 and I2.4.5[ 

2.4.3.1 Virtual corrections: UV singularities and counterterms 

Self-energy and vertex one-loop corrections to the tree level gg — * Zbb process (see 
Figures 12.161 and I2.17P give rise to UV divergences. These singularities are canceled by 
a set of counterterms fixed by well defined renormalization conditions. We need to introduce 
counterterms for the external field wave functions of bottom quarks and gluons {5Z2^\ SZ^), 
for the bottom mass {Snih), and for the strong coupling constant {dZ^J. 

By carefully grouping subsets of self-energy and vertex diagrams, we can factor out the 
UV singularities of the 0{al) virtual amplitude and write them in terms of the tree level 
partial amplitudes Aq^s, -^0,1, and Aq^u introduced in Eq. (12.561) . According to the notation 
introduced in Figures 12. 16112. 11J| we denote by Dij (with D = S,V, i = l,2, . . ., and j = s, t, u) 
a class of diagrams with a given self-energy or vertex correction insertion, summed over all 
possible insertions of the external Z boson, one for each different diagram. We now define 
Auv(^Dij ) to be the UV pole part of the corresponding amplitude. Using this notation, we 
find: 
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(2.67) 



where nif = 5 corresponds to the number of hght quark flavors, A^ = 3 is the number of colors 
and Ms and Afb are defined in Eq. (12.191) . 

We notice that some of the UV divergent virtual corrections (V^i '^i]^i^u) ^ ^^"^ ^''(fu")^)' 
as well as 6Z2^^ and 6Z3 in Eqs. (12.201) and (I2.68P below, have also IR singularities. In this 
section, we limit the discussion to the UV singularities only, while the IR structure of these 
terms will be considered in Section I2.4.3.2[ For this reason, we have explicitly denoted by 
euv the pole parameter. 

The corresponding counterterms are defined as follows. For the external fields, we fix the 
wave-function renormalization constant of the external bottom quark fields using the on-shell 
subtraction scheme, giving the expressions shown in Eq. (12.201) . while we renormalize the 



wave-function of external gluons in the MS subtraction scheme: 
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according to which we also need to consider the insertion of a finite self-energy correction on 
the external gluon legs. This amounts to an extra contribution 
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(2.69) 



which is important in order to obtain the correct scale dependence of the NLO cross section. 

The structure of the remaining counterterms that appear in this case is the same as 
presented in Sections 12.3.1.11 and 12.4.2. 1^ and we refer the reader to them. 

Using the results in Sections 12.3.1.11 and I2.4.2.T] and in Eqs. (I2.67p - fl2.69p . it is easy to 
verify that the UV pole part of cr^g^g^, 



(2.70) 



is free of UV singularities and has a residual renormalization scale dependence of the form: 

as expected by renormalization group arguments (see the first term of Eq. (12. 7p ). We note 
that the presence of s in the argument of the logarithm of Eq. fl2.7ip has no particular 
relevance as described at the end of Section [2.3.1.11 

2.4.3.2 Virtual corrections: IR singularities 

The structure of the IR singularities originating from the 0{as) virtual corrections to the tree 
level amplitude for gg Zbb is more involved than that for the UV singularities. However 
it simplifies considerably when given at the level of the amplitude squared, and this is what 
we present in this section. 

The IR divergent part of the 0{al) virtual amplitude squared of Eq. (12.651) can be written 
in the following compact form: 



J2 E^^^ ^-(^A.)) = ^^^E (^i-^S + + CsM^vl) , (2.72) 



where Aft is defined in Eq. (12.191) and we denote by Air(^£i^^) the IR pole part of the 
amplitude of a given Di j class of diagrams. The result is organized in terms of leading and 
sub- leading color factors presented in Eq. (I2.66p . and the corresponding matrix elements 
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squared My], My], and Myj^ are given by: 
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where the IR nature of the pole terms has been made exphcit. A'^ and Aq""^ are defined in 
Eq. ( 12.56p . while v4.o,s, ^o,*, and Aq^u are given explicitly in Appendix [Bl We have introduced 
the notation Ag = ln(s/m^) and A^-^ = ln(rj/m^), where the invariants have been defined in 
Eq. ([22SD, and (3^,i and A^j are defined in Eq. flZ^ . 

We write explicitly here the IR divergent structure of the {5Z-i) counterterm, as it is the 
only one that has not appeared previously. We get: 



^2.74) 



When we add the IR singularities coming from the counterterms we have introduced in 
Section [2.4.3. ![ we can write the complete pole part of the IR singular 0{al) virtual cross 
section as: 



( ;^virt \ 



gg )IR-pole 



a. 



1 \ 1 



(2.75) 



27r" V3 3" ' N/ CiR, 
As will be demonstrated in Section 12.4.41 the IR singularities of dgg^ are canceled by the 
corresponding IR singularities of d'g^^^. 

2.4.4 Real Corrections to Zbb Production 

The real corrections to Zbb production is built up of three 0{al) subprocesses, namely 
qq Zbb + g, gg Zbb + g and (g, q)g Zbb + (g, q). 



53 




Figure 2.20: Examples of 0{as) real corrections to qq Zbb production, with Z emitted 
from the final fermion line. 
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Figure 2.21: Examples of 0{as) real corrections to Zbb production. The circled crosses 
denote all possible insertions of a Z weak boson, each insertion corresponding to a different 
diagram. 



The real cross section for the qq Zbb + g subprocess is analogous to the qq' — > Wbb + g 
subprocess, although it contains a set of diagrams with the Z boson emitted from the 6-quark 
fermion line, as shown in Figure I2.20[ We refer the reader to the discussion in Sections 12.3.21 
and 12.3.31 for details. 

The NLO real cross section a^^^^^ in Eq. (12.51) corresponds to the 0{as) corrections to 
gg Zbb due to the emission of a real gluon, i.e. to the process gg Zbb + g, examples of 
which are illustrated in Figure 12.211 It contains IR singularities that cancel the analogous 
singularities present in the 0{as) virtual corrections (see Section 12.4.3.21) and in the NLO 
parton distribution functions. These singularities can be either soft, when the energy of 
the emitted gluon becomes very small, or collinear, when the final state gluon is emitted 
coUinear to one of the initial gluons. There is no collinear radiation from the final b and b 
quarks because they are massive. At the same order in a^, the 0"^°^' cross section corresponds 
to the tree level processes (g, q)g Zbb + (g, q), an example of which is also illustrated in 
Figure I2.21[ This part of the NLO cross section develops IR singularities entirely due to the 
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coUinear emission of a final state quark or antiquark from one of the initial state massless 
partons. The IR singularities can be conveniently isolated by slicing the gg — ^ Zbb + g and 

'l)g ~^ Zbb + (g, q) phase spaces into different regions defined by suitable cutoffs (as was 
done in Section r2.3.2p . The dependence on the arbitrary cutoffs introduced in slicing the 
phase space of the final state particles is not physical, and cancels at the level of the total 
real hadronic cross section, i.e. in a'^'^^^, as well as at the level of the real cross section for each 
separate channel, i.e., in cr^|'^\ ^IT^^ ^^"^ '-'^IT^' This cancellation constitutes an important 
check of the calculation and will be discussed in detail in Section I2.4.5[ 

We have calculated the cross section for the processes 

9{qi) + 9{q2) b{pb) + b{pi) + Z{pz) + g{k) , 

and 

{q, + 9(^2) b{pb) + b{pi) + Z{pz) + (g, q){k) , 

with qi + q2 = Pb + Vb + Vz + k, using the PSS method with two cutoffs, as in Section 12.3.21 
In the next section we discuss the application of the two-cutoff PSS method to the 
gg Zbb + g and q{q)g Zbb + q{q) subprocesses. 

Real gluon emission, gg Zbb + g: soft region 

The soft region of the phase space for the gluon emission process 

9\qi) + 9\q2) ^ b{p,) + b{pi) + Z{pz) + /(A;) (2.76) 
is defined by demanding that the energy of the emitted gluon {k^ = Eg) satisfies the condition 

Eg < (2.77) 

for an arbitrary small value of the soft cutoff 5s- In the soft limit {Eg — * 0), the amplitude 
for this process can be written as: 

Asofti.99 Zbb + g) = 

- (1^ - 1^ ) <-^- ^ - (1^ - IS) <-^- - 

+ tvf (*:£ - *±) (A, + A,,.) + (^ - ^) (A,,. - A,.) . 

\qi-k q2-k J \q2-k qi-k J 

(2.78) 
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where a, b, and c are the color indices of the external gluons, while e'^{k, A) (for A = 1, 2) is the 
polarization vector of the emitted soft gluon. Moreover, in the soft region the gg Zbb + g 
phase space factorizes as: 

d{PS,)igg^Zbb + g) ^ d{PSs){gg ^ Zbb)d{PSg)soft 



(277) ('^-1)2^/ V 2 



d{PSs){gg - Zbb) ^ 9 ( 6,^ - E, 



(2.79) 



where d{PS4) and d^PS^) have been defined in Section 12.11 while d{PSg)soft denotes the 
the phase space measure of the soft gluon. Since the contribution of the soft gluon is now 
completely factorized, we can perform the integration over d{PSg)soft analytically, using 
dimensional regularization in d = 4 — 2e to extract the soft poles that will cancel the 
corresponding singularities in Eq. (12.751) . For completeness, the integrals that we have used 
to perform the integration over the phase space of the soft gluon are collected in Appendix [Fl 
After squaring the soft amplitude Asoft, summing over the polarization of the radiated 
soft gluon, and integrating over the soft gluon momentum, the pole part of the parton level 
soft cross section reads 



pole 



(^^f )po/e = J d{PS,) (^j d{PSg)softY.\Asoft{gg ^ Zbb + g) 

= j d{PS,)^U^{c,M^l + C,M^l + C,Mf) , (2.80) 

where Ci, C2, and C3 are defined in Eq. (12.661) . while M.'s\-, ■M.'s\i ^"^^ -^s^e represent the 
IR pole parts of the corresponding matrix elements squared, and can be written as: 

2 



-^S^ = -M'^il--^{l+AH5s)){\AT''? + \Af\'') 



Mf^ = -Mf, + -\Af\^. (2.81) 



e 



Note that in this section we do not explicitly denote the IR poles as poles in eia, since all 
singularities present in cr^gg^g are of IR origin. 

After adding the IR divergent part of the parton level virtual cross section of Eq. (12.751) 
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we obtain: 



99 



soft- 
99 ' 



pole 



V^gg JIR^pole 



-4iV In (5, 



-(UN 
3^ 



2n 



If) 



—a 



LO 
99 



(2.82) 



where we can see that the IR poles of the parton level virtual cross section are exactly canceled 
by the corresponding singularities in the parton level soft gluon emission cross section. The 
residual divergences will be canceled by the soft+virtual part of the PDF counterterm when 
convoluting with the gluon PDFs, as will be demonstrated in Section 12.4.51 The finite 
contribution to the parton level soft cross section is finally given by 



^lf)fimte = / diPSs) I / d{PSg)softY,\'^soft{g9 ^ Zbb + g)\' 



finite 



(2.83) 



where the finite parts of the , and Ai^^' matrix element squared are explicitly 

given by: 



If. 



,(3) 
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(2) 
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(3) 



where 



—n^ + iA, ln(5,) +81n2((5,) - 2A, -41n((5,) 
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1-A 
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1 1 

(A, + 21n(5,)) (A., + + -F{qi,pb) + -F{q2,Pi) 
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^P'bil-P^) 
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p-b{^-^) 
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2A, - 41n(5,) + {\AT'? + 

-^TT^ - 2A2 + 8 ln2(5,) + 2 (A, + 2 ln(5,)) (A,, + A,, + A,3 + A, 
+ F{q2,pi) + Ffe,^;,) + F{q^,pi) 
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1 + A: 



1) 



(2.84) 



(2.85) 



(2.86) 



l-A-6 ' "^"^ 2(a,-,pO-p5) ' 
A^ = ln((5/m^), with 5 a given kinematic invariant, and 5^,5, AE; A^j are defined in 
Eq. (12.261) . while the function F{pi,pf) can be found in Eq. (IF.Op (see Appendix [Fl for more 
details) . We have used the set of kinematic invariants presented in Eq. (12.231) . 
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Real gluon emission, gg — > Zhh + g: hard region 

The matrix element squared for ij Zbb + g factorizes into the LO matrix element squared 
and the unregulated Altarelli-Parisi splitting function Pa/ = P^-, + eP/-, for i — > i'g, i.e.: 



ElAcai(^J - Zbb + g)\' ''"'^ (4vra.)5^|^o(^'j - Zbb)\' 



z s 



(2.87) 



«3 



where Sig = 2qi-k and P^^, and P/j, are the Altarelli-Parisi splitting functions, which in the gg 



case are given by: 



2N 



0. 



1 - z 



1 - ^ 



z(l 



{21 



Moreover, in the collinear limit, the ij —>■ Zbb + g phase space also factorizes as: 
d{PS,){zj ^ Zbb + g) d{PS,){i'j Zbb) e (e„- 5,,^ I X 



^(cose,,-(l-4)) 



d=4-2e 



d{PS^){i'3^Zbb) 



— \ Z dz dSia f(l — z)Sia] 

r 1 - e 167r2 ^ 



X 



(2.89) 



where d{PS4) and d{PSs) have been defined in Sect ion [2Tn while the integration range for Sig 
in the collinear region is given in terms of the collinear cutoff, and we have defined s' = 2qii-qj. 
The integral over the collinear gluon degrees of freedom can then be performed analytically, 
and this allows us to extract explicitly the collinear singularities of d'gg''''^ [HI [91], which can 
be written as: 

hard/ coll 



99 



a. 



4TTfl 

27rr(l -e) V^«r 



2 \ ^■ 



dz 



[l-zf s' 
2z ml 



67' X 



Pu'{z)a^;;{z'j^Zbb) + {z^j) 



[2.90) 



where , and j are all gluons. As usual, these initial state collinear divergences are absorbed 
into the gluon distribution functions as will be described in detail in Section 12.4.51 
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The tree level processes (g, q)g — > Zbb + {q, q) 

The extraction of the collinear singularities of 0"^°^' is done in the same way as described 
in the previous subsection for the gg initial state. In the collinear region, cos > 1 — 5ci 
the initial state parton i with momentum qi is considered to split into a hard parton i' and 
a collinear quark q, i i'q, with g^/ = zqt and k = {1 — z)qi. The matrix element squared 
for ij — s> Zhh + q factorizes into the unregulated Altarelli-Parisi splitting functions in d 
dimensions, Pw = P^^, + eP-^,, shown in Eq. (12.471) for the g ^ qq splitting and in Eq. (12.921) 
for the {q,q) g{(l,(j) splitting, and the corresponding LO matrix elements squared. The 
ij Zbb + q phase space factorizes into the i'j Zbb phase space and the phase space of 
the collinear quark. As a result, after integrating over the phase space of the collinear quark, 
the collinear singularity of a^^^^j can be extracted as: 

,2^ 



±.coll 



as 1 /47r^^^^^ ' 



67' / dz 



z? s' 



2z ml 



X 



_27rr(l -e) \ ml 

[P^oiz) al^{9{qi')9{q2) - Zbb) + P,,{z) a^^^iqiq^Mq,,) ^ Zbb)] . (2.91) 

The collinear radiation of an antiquark in qg —* Zbb + g is treated analogously. The 0(1) 
and 0(e) parts of the Pgg splitting function are given by: 

^;(^) = -Cfz, (2.92) 

and by Eq. (12:471) for Pgg. 

Again, these initial state collinear divergences are absorbed into the parton distribution 
functions as will be described in detail in Section 12.4.51 

2.4.5 Total Cross Section of pp{pp) Zbb 

The total inclusive hadronic cross section for pp{pp) — * Zbb is the sum of the contributions 
from the qq, the gg and the (g, q)g initial states: 

^ Zbb) = a^,^°{pp{pp) ^ Zbb) + a-g^'^ippipp) ^ Zbb) + C°(pp(pp) ^ Zbb) . 

(2.93) 

As described in Section I^TTl crfj^^ippipp) Zbb) is obtained by convoluting the parton level 
NLO cross section d-fj^° {pp{pp) Zbb) with the NLO PDFs J-'f'''{x, fi) {i = q,g), thereby 
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absorbing the remaining initial state singularities of SafJ"'^ into the renormalized PDFs. In 
the following we demonstrate in detail how this cancellation works in the case of the gg and 
(g, q)g initiated processes. The case of the qq initiated process is discussed in Section 12.3. 3^ 
where we presented in detail the contribution of the qq initial state to a^^°{pp{pp) — Wbb). 
'^q^'^ {pp{pp) ~^ Zhh) can be obtained from there with obvious modifications, and will not be 
repeated here. 

Similarly to the discussion in Section [2.3.3[ first the parton level cross section is convoluted 
with the hare parton distribution functions '^(x) and subsequently the J-'f'^{x) are replaced 
by the renormalized parton distribution functions, J^f'^{x, fif), defined in some subtraction 



scheme at a given factorization scale /ij. Using the MS scheme, the scale-dependent NLO 
parton distribution functions for the two-cutoff PSS are given in terms of the bare JF^ 



■p,Pi 



X] 



and the QCD NLO parton distribution function counterterms [79] as follows: 
(a) For the case where an initial state gluon, quark or antiquark {k = g,{q,q)) splits. 



respectively, into a qq or (g, q)g pair {k' = (g, q),g): 



27r \ I r(l-e] 



^ dz 



■p,p 



X 



(2.94) 



where P^j are defined in Eqs. (I2.47p and (I2.92p 
(b) For the case of g ^ gg splitting: 



T^'^ix) 



271 I jjj 




ri 



1-Ss 



(2.95) 



where Pgg is Altarelli-Parisi splitting function given in Eq. fl2.88p . 



The 0{as) terms in the previous equations are calculated from the 0{as) corrections 
to the k k'j splittings, in the two-cutoffs PSS formalism. Moreover, note that in 
Eqs. (12.941) and (I2.95P we have carefully separated the dependence on the factorization (///) 
and renormalization scale (/i^). It is understood that = as{fir)- The definition of the 
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subtracted PDFs is indeed the only place where both scales play a role, and the only place 
where we have a dependence on ^f. In the rest of this section, we have always set 
and we will also give the master formulas for the total NLO cross section, Eq. f l2.96p . using 
In Chapter [3] we will study the explicit dependence of the total hadronic cross 
section on both scales /i^ and /i/. 

When convoluting the parton gg cross section with the renormalized gluon distribution 
function of Eq. fl2.95p . the IR singular counterterm of Eq. fl2.95l) exactly cancels the remaining 
IR poles of a^g]^^ + a'f/* and a-^^'^/^"", shown in Eq. Finally, the complete O(a^) 

inclusive total cross section for pp{pp) Zbb in the MS factorization scheme when only the 
gg initial state is included, i.e. cr^q^ipvivv) ~^ Zbb) of Eq. fl2.93p . can be written as follows: 



NLO 
99 



= \ j dXldX2 {:F^{xi, fl)J^^P\x2, fl) [^gg{xi, X2, fl) + {alf)f^nite{xi,X2, /i) 
+ {^gg^')f^nite{Xl, X2, fi) + <J^+^'+^* (Xl , X2 , /i) + (1 ^ 2)] } 

X a-(x„x„/.)^ln (^^^^1) P99i^) + (1 - 2)} 



2)} , 



(2.96) 



where <Tg+''+^* is obtained from the sum of (<3-J^*)c/y-poZe of Eq. fl2T0D . of Eq. I^M), 



and the PDF counterterm in Eq. fl2.95p as follows 



a 



s+v+ct 
99 
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2^ 



4iVln((5,)ln 



1^' 
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'-N 



In. 



If 



g +4iVln(5,)jln( ^ 
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99 



(2.97) 



We note that 0"^^° is finite, since, after mass factorization, both soft and coUinear 



singularities have been canceled between a™^ + a^"^^ and a. 
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hard/ coll 
'99 



The last terms 



respectively describe the finite real gluon emission of Eq. 02.121) . Note that when collecting 
all the terms in Eq. 02.960 that are proportional to ln(/i^/s), one obtains exactly the last two 
terms in Eq. 02. 7p . as predicted by renormalization group arguments. 
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For the (g, q)g initiated processes we find 



19 
NLO 



99 ' 



Xi,X2,/i) 
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1 dz 
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^f(-,/x)^f^Hx2,/i)x 
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LO, 
99 



p;,(^)in 
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i=q,q 



+ (1^2) 



2)} 



(2.98) 



We would like to conclude this section by showing explicitly that the total NLO cross 
section, a^^°, does not depend on the arbitrary cutoffs introduced by the PSS method, i.e. 
on 5s and 5c- The cancellation of the PSS cutoff dependence is realized in a'''^'^' by matching 
contributions that are calculated either analytically, in the IR-unsafe region below the cutoffs, 
or numerically, in the IR-safe region above the cutoffs. While the analytical calculation in 
the IR-unsafe region reproduces the form of the cross section in the soft or collinear limits 
and is therefore only accurate for small values of the cutoffs, the numerical integration in the 
IR-safe region becomes unstable for very small values of the cutoffs. Therefore, obtaining a 
convincing cutoff independence involves a delicate balance between the previous antagonistic 
requirements and ultimately dictates the choice of values that are neither too large nor too 
small for the cutoffs. The Monte Carlo phase space integration has been performed using 
the adaptive mult i- dimensional integration routine VEGAS [7T] . 

Figures [2.221 and 12.231 illustrate the dependence of the total cross section o"'^^°(pp Zbb) 
on the two-cutoffs of the PSS method, using the setup outlined in Section [XTl In Figure U.2'2\ 
5s is varied between 10"^ and 10"^ with 5^ = 10"^, while in Figure [2.231 5c is varied between 
10"'' and 10~^ with 5s = 10~^. In both plots, we show in the upper window the overall 
cutoff dependence cancellation between a'^"-^^ + a'^^rd/coa ^^^^ ^hard/non-coii ^reai including 

all channels, gg, qq and qg. We include too contributions from the LO and the virtual cross 
sections which are cutoff independent. In the lower window of the same plots we complement 
this information by reproducing the full cr^^°, including all channels, on a larger scale that 
magnifies the details of the cutoff dependence cancellation. The statistical errors from the 
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Figure 2.22: Dependence of a'^^°{pp Zbb) on the soft cutoff 6s of the two-cutoff PSS 
method for fi = Mz + 2mb, and 6c = 10"^. The upper plot shows the cancellation of the 6s- 
dependence between Q-^oft _^^hard/coU ^hard/non-coii ^ rpj^^ lower plot shows, on an enlarged 

scale, the dependence of the full a^^'^ 
statistical errors. 



+ 0"^^° + ° on 6s with the corresponding 
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Monte Carlo phase space integration are also shown. Both Figures 12.221 and 12.231 show a 
clear plateau over a wide range of 6s and 6c and the NLO cross section is proven to be cutoff 
independent. The results presented in Chapter [3] have been obtained by using (5s = 10~^ and 
6c = 10-'. 
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Figure 2.23: Dependence of a^^°{pp —* Zbb) on the coUinear cutoff 6c of the two-cutoff 
PSS method, for /i = + 2mb, and 5, = IQ-^ The upper plot shows the cancellation of 
the 5s-dependence between a^°-^^ + o-^'^'^'^/co"^ and (jhard/non-coii ^ r^^^Q lower plot shows, on 



an enlarged scale, the dependence of the full (T^^° 
corresponding statistical errors. 
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CHAPTER 3 



Numerical Results 



In this Chapter, we present numerical results for the total cross sections and distributions 
for W/Z bb production including NLO QCD corrections and complete bottom-quark mass 
effects [m [92]. We specialize our discussion to the case of the Tevatron collider, because 
this is at the moment the most interesting phenomenological environment (see introduction 
in Chapter [1]). We also investigate the stability of the NLO QCD results by studying the 
dependence of the total cross section on the renormalization [fir) and factorization (fif) 
scales. Finally, we carefully compare our results with results obtained from a NLO calculation 
that considers massless bottom quarks by using the MCFM code ^9] . 

In Section 13.11 we specify the setup used to produce the plots, while in Sections 13.21 and 
13.31 we present and discuss results for Wbb and Zbb production respectively. 

3.1 The Setup 

We present NLO QCD results for W/Z bb production at the Tevatron using a non-zero 
bottom-quark mass fixed at mb=4.62 GeV. The W and Z bosons are considered on-shell 
and their masses are taken to be Mw = 80.410 GeV and Mz = 91.1876 GeV. The 
mass of the top quark, entering in virtual corrections, is set to nit = 170.9 GeV. The 
LO results use the one-loop evolution of as and the CTEQ6L1 set of PDF [93], with 
a^'^{Mz) = 0.130 , while the NLO results use the two-loop evolution of as and the CTEQ6M 
set of PDF, with a^^'~'{Mz) = 0.118. The W boson coupling to quarks is proportional to 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. We take Vud = = 0.975 and 
Vus = Vcd = 0.222, while we neglect the contribution of the third generation of quarks, since 
it is suppressed either by the initial state quark densities or by the corresponding CKM 
matrix elements. 
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Table 3.1: LO and NLO total Wbb cross sections at the Tevatron for massive and massless 
bottom quarks, using fir = f^f = Mw + 2m;,. The numbers in square brackets are the ratios 
of the NLO and LO cross sections, the so called A'-factors. Statistical errors of the MC 
integration amount to about 0.1%. 



Cross Section 


mi, (pb) [ratio] 


nib = (pb) [ratio] 


0"LO 


2.20[-] 


2.38[-] 


o"nlo inclusive 


3.20[1.45] 


3.45[1.45] 


o"NLO exclusive 


2.64[1.20] 


2.84[1.19] 



Partons cannot be detected as they are always confined in hadrons. For this reason, 
any phenomenological collider study, including final hadronic states, must implement a 
jet algorithm to recombine partons into jets, in a way consistent with factorization (see 
Section 12.11) and with experimental techniques. The jet algorithm basically assigns a 
"separation" between partons, and, based on it, defines criteria to decide whether to group 
a set of partons in the final state into a "proto-jet". Finally, kinematic cuts are applied 
depending on the experimental setup, to decide if a proto-jet is in an observable region, in 
which case the proto-jet is promoted to a jet. We will consider in our study 6-type and 
light-type jets, where the former contains either a 6 or a 6 quark and the latter can only 
contain massless quarks or gluons. 

We implement the fcr jet algorithm [94l [95l [96l [97] with a pseudo-cone size R = 0.7 and we 
recombine the parton momenta within a jet using the so called covariant i?-scheme [95]. We 
checked that our implementation of the k-r jet algorithm coincides with the one in MCFM. 
We require all events to have a bb jet pair in the final state, with each jet having a transverse 
momentum larger than 15 GeV (p^** > 15 GeV) and a pseudorapidity that satisfies \r]^'^\ < 2. 
We impose the same px and |?7| cuts also on the extra jet that may arise due to hard non- 
coUinear real emission of a parton, i.e. in the processes W/Z bb + g or W/Z bb + q{q). This 
hard non-coUinear extra parton is treated either inclusively or exclusively. In the inclusive 
case we include both two- and three-jet events, while in the exclusive case we require exactly 
two jets in the event. Two-jet events consist of a bottom-quark jet pair that may also include 
a final-state light parton (gluon or quark) due to the applied jet algorithm. Results in the 
massless bottom-quark approximation have been obtained using the MCFM code |49j . 
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Figure 3.1: Dependence of the LO (black solid band), NLO inclusive (blue dashed 
band), and NLO exclusive (red dotted band) Whh total cross sections on the renormaliza- 
tion/factorization scales, including full bottom-quark mass effects. The bands are obtained 
by independently varying both //^ and /i/ between /io/2 and 4/io (with = m}y + Mw /2). 



3.2 Whh Production at the Tevatron 

Let us first consider the influence of the NLO QCD corrections on the total cross section. 
In Table 13.11 we present the results obtained for both LO and NLO total cross sections, at 
a reference scale fir = fJ^j = My/ + 2m;,, both in our fully massive calculation and in the 
massless approximation. 

It can be seen that, given the setup explained in Section [3?T| the NLO QCD corrections 
increase considerably the total cross section, with NLO/LO ratios (/^-factors) of about 1.45 
and 1.2 for the inclusive and exclusive case respectively (for both the massive and massless 
calculations). We can see also that, in general, the massless approximation overestimates 
the total cross section. In the following we will study in detail where these corrections are 
more important, and we will show that, in the case of distributions, a global rescaling (or 
i^-factor) does not properly simulate the NLO corrections. 

In Figures 13.11 and 13.21 we illustrate the renormalization and factorization scale depen- 
dence of the LO and NLO total cross sections, both in the inclusive and exclusive case. 
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Figure 3.2: Dependence of the LO and NLO inclusive and exclusive pp Wbb total cross 
section on the renormalization/factorization scale, when fir = fJ'f = fi- The left hand side 
plot compares both LO and NLO total cross sections for the case in which the bottom 
quark is treated as massless (MCFM) or massive (our calculation). The right hand side plot 
shows separately, for the massive case only, the scale dependence of the qq' and qg + qg 
contributions, as well as their sum. 
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Figure 3.3: Dependence on the renormalization/factorization scale of the rescaled difference 
between our NLO calculation (with rrib ^ 0) of the total Wbb cross section and the 
corresponding result computed using MCFM (with = 0) for the inclusive and exclusive 
cases (with /ir = /i/ = /io) respectively. The error bars indicate the statistical uncertainty of 
the Monte Carlo integration. 



Figure [XT] shows the overall scale dependence of both LO, NLO inclusive and NLO exclusive 
total cross sections, when both fij. and fif are varied independently between /io/2 and 4/io 
(with fiQ = rrih + M\y/2), including full bottom-quark mass effects. We notice that the 
NLO cross sections have a reduced scale dependence over the range of scales shown, and the 
exclusive NLO cross section is more stable than the inclusive one especially at low scales. 
This is consistent with the fact that the inclusive NLO cross section integrates over the entire 
phase space of the qg{qg) bbW + q{q) channels that are evaluated with NLO and NLO 
PDFs, but are actually tree-level processes and retain therefore a strong scale dependence. 
In the exclusive case only the 2 — 3 collinear kinematic of these processes is retained, since 
3-jets events are discarded, and this makes the overall renormalization and factorization scale 
dependence milder. To better illustrate this point, we show in the right hand side plots of 



Figures 3.2(a) and 3.2(b) the /i-dependence of the total cross section and of the partial cross 
sections corresponding to the qq' and the qg + qg initiated channels separately, for fir = fif, 
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both for the inclusive and for the exclusive case. It is clear that the low scale behavior of 
the inclusive cross section is mainly driven by the qg + qg contribution. In the left hand 



side plots of Figures 3.2(a) and 3.2(b), we also compare the scale dependence of our results 
to the scale dependence of the corresponding results obtained with m^, = (using MCFM), 
both at LO and at NLO. Using a non-zero value of mh is not expected to have any impact 
on the scale dependence of the result^] and, indeed, the scale dependence of the LO and NLO 
pair of curves is very similar, with a shift due to the bottom-quark mass effects. 

While the LO cross section still has a 40% uncertainty due to scale dependence, this 
uncertainty is reduced at NLO to about 20% for the inclusive and to about 10% for the 
exclusive cross sections. The uncertainties have been estimated as the positive/negative 
deviation with respect to the mid-point of the bands plotted in Figure 13.11 where each 
band range is defined by the minimum and maximum value in the band. We notice 
incidentally that the difference due to finite bottom-quark mass effects is less significant 
than the theoretical uncertainty due to the residual scale dependence in the inclusive case, 
but is comparable in size in the exclusive case. Indeed, the finite bottom-quark mass effects 
amount to about 8% in both inclusive and exclusive cases. 

In Figure 13.31 we show the rescaled difference between the NLO total cross sections 
obtained from our calculation (with mf, ^ 0) and with MCFM (with m;, = 0) defined as 
follows: 

Aa = + 0) - = 0) "^^^^^ ■ (3.1) 

As can be seen, within the statistical errors of the Monte Carlo integration, the finite bottom- 
quark mass effects on the total cross sections at NLO are well described by the corresponding 
effects at LO. 

Finally, in Figures [3^3.61 we study the distribution dajdm^i^ where m^j is the invariant 
mass of the 66 jet pair. The impact of NLO QCD corrections on this distribution is illustrated 



in Figures 3.4(a) and 3.4(b) for the inclusive and exclusive case, respectively. We see that the 
NLO QCD corrections affects the cross section quite substantially in particular for low values 
of mfog. In each figure the right hand side plot gives the ratio of the NLO and LO distributions. 



providing a sort of K-factor bin by bin. Figures 3.5(a) and 3.5(b) compare the NLO da jdm^^i 
distributions obtained from the massive and massless bottom-quark calculations. The results 



^Note that we always use mh = 4.62 GeV in the determination of the scales in terms of /io = m-6 + 
even in the results obtained with rub = 0. 
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i^bb (GeV) m^, (GeV) 

(b) Exclusive case 

Figure 3.4: The distribution da{pp Wbb)/dmf^i in LO and NLO QCD. The right hand 
side plot shows the ratio of the LO and NLO distributions. 
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(b) Exclusive case 



Figure 3.5: The inclusive and exclusive distributions da{pp — ^ Whb) / dm^i derived from our 
calculation (with mi, ^ 0) and from MCFM (with = 0). The right hand side plot shows 
the ratio of the two distributions, da{mb 7^ 0)/da{mh = 0). 
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Figure 3.6: The LO distribution da{pp Wbb)/dmiji derived from our calculation (with 
rrih 7^ 0) and from MCFM (with mi, = 0). The right hand side plot shows the ratio of the 
two distributions, da^niiy ^ 0)/da{'mb = 0). 



pp Wbb 




m,j(GeV) 



Figure 3.7: The mf,i distribution of the rescaled difference between our NLO calculation 
(with nih 7^ 0) and MCFM (with nib = 0) for the inclusive (upper plot) and exclusive (lower 
plot) pp Wbb production. 
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with rrify = have been obtained using MCFM. As expected, most of the difference between 
the massless and massive bottom-quark cross sections is coming from the region of low 
invariant mass m^^;, both for the inclusive and exclusive case, where the cross sections for 
7^ are consistently below the ones with = 0. For completeness, we also show in 
Figure 13.61 the comparison between massive (m^ 7^ 0) and massless (m^ = 0) calculations 
at LO in QCD. The LO rrif^i distribution for massive bottom-quarks has been obtained 
both from our calculation and from MCFM, which implements the m^, ^ option at tree 
level, and both results have been found in perfect agreement. As can be seen by comparing 



Figures 3.5(a) 3.5(b) and Figure ESI the impact of a non-zero bottom-quark mass is almost 
not affected by including NLO QCD corrections. To illustrate this in more detail, we show 
in Figure 13.71 the rescaled difference between the m^g distributions obtained with our NLO 
calculation (with 7^ 0) and with MCFM (with rrii, = 0) defined as follows: 

We notice that finite bottom-quark mass effects are particularly relevant for values of the 
invariant mass below about 60 GeV and that they appear to be of the same order at 
LO and NLO. 

3.3 Zbb Production at the Tevatron 

To start, let us have a look at the influence of the NLO QCD corrections on the total cross 
section. In Table [312] we present the values obtained with the scale fir = f^f = Mz + 2mh, 
considering LO and NLO total cross sections, both in our fully massive calculation and in 
the massless approximation. 



It can be seen that, given the setup explained in Section \3A] the NLO QCD corrections 
increase considerably the total cross section, with NLO vs. LO ratios (i^'-factors) of about 1.5 
and 1.27 for the inclusive and exclusive case respectively (for both the massive and massless 
calculations). We can also see that, in general, the massless approximation overestimates 
the total cross section. In the following we will study in detail where these corrections are 
more important, and especially we show that, in the case of distributions, a global rescaling 
(or A'-factor) does not properly simulate the NLO corrections. 

In Figures ESI and ES] we illustrate the renormalization and factorization scale dependence 
of the LO and NLO total cross sections, both in the inclusive and exclusive case. Figure 13.81 
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Figure 3.8: Dependence of the LO (black solid band), NLO inclusive (blue dashed 
band), and NLO exclusive (red dotted band) Zhh total cross sections on the renormaliza- 
tion/factorization scales, including full bottom-quark mass effects. The bands are obtained 
by independently varying both /i.^ and yU/ between /io/2 and 4/io (with fiQ = nih + Mz/2). 



shows the overall scale dependence of both LO, NLO inclusive and NLO exclusive total 
cross sections, when both /x^ and fif are varied independently between /io/2 and 4/io (with 
fiQ = rrih + Mz/2), including full bottom-quark mass effects. We notice that the NLO cross 
sections have a reduced scale dependence over the range of scales shown, and the exclusive 
NLO cross section is more stable than the inclusive. Similarly to what we have discussed in 
the Wbb case, this effect is mainly driven by the tree level subprocess q{q)g Zbb + q{q) 
contributing to the real corrections. In the Zbb case, we also have a new initial state, namely 
gg. Its scale dependence behavior is similar to the qq initiated subprocess. To illustrate the 



independent contributions, we show in the right hand side plots of Figures 3.9(a) and 3.9(b) 



the /i-dependence of the total cross section and of the partial cross sections corresponding to 
the qq, qg + qg and gg initiated channels separately, for fir = fJ'f, both for the inclusive and 
for the exclusive case. It is clear that the low scale behavior of the inclusive cross section 
is considerably affected by the qg + qg contribution, which show a monotonic dependence 
on /i (i.e. with no plateau) characteristic of tree level processes. In the left hand side 



plots of Figures 3.9(a) and 3.9(b) we also compare the scale dependence of our results to 
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Table 3.2: LO and NLO total Zbb cross sections at the Tevatron for massive and massless 
bottom quarks, using fir = fJ'f = Mz + 2mb. The numbers in square brackets are the ratios 
of the NLO and LO cross sections, the so called i^-factors. Statistical errors of the MC 
integration amount to about 0.1%. 



the scale dependence of the corresponding results obtained with = (using MCFM), 
both at LO and at NLO. Using a non-zero value of mb is expected to have a mild impact 
on the scale dependence of the results, as the only modification to the renormalization scale 
dependence comes from the bottom quark mass renormalization, as shown in the subsections 
of Section [23] in Chapter [2l Indeed, the scale dependence of the LO and NLO curves is very 
similar. 

While the LO cross section still has a 45% uncertainty due to scale dependence, this 
uncertainty is reduced at NLO to about 20% for the inclusive and to about 11% for 
the exclusive cross sections. As before, the uncertainties have been estimated as the 
positive/negative deviation with respect to the mid-point of the bands plotted in Figure 
where each band range is defined by the minimum and maximum value in the band. 
We notice incidentally that the difference due to finite bottom-quark mass effects is less 
significant than the theoretical uncertainty due to the residual scale dependence in the 
inclusive case, but is comparable in size in the exclusive case. Indeed, the finite bottom- 
quark mass effects amount to a reduction of the total cross sections by about 7% compared 
to the massless case at both LO and NLO QCD. 

In Figure [3.101 we show the rescaled difference between the total cross sections obtained 
from our calculation (with rrib 7^ 0) and with MCFM (with rrib = 0) defined as in Eq. (13. ip . 
As can be seen, within the statistical errors of the MC integration, the finite bottom-quark 
mass effects on the total cross sections at NLO are well described by the corresponding effects 
at LO, similarly to what is observed in the Wbb case. 

Finally, in Figures 13.111 to 13.131 we study the distribution dcr/dnibi, where rUbb is the 
invariant mass of the bb jet pair. The impact of NLO QCD corrections on this distribution 



Cross Section 
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Figure 3.9: Dependence of the LO and NLO inclusive and exclusive pp Zbb total cross 
section on the renormahzation/factorization scale, when fir = fJ'j = fJ'- The LHS plots 
compare both LO and NLO total cross sections for the case in which the bottom quark is 
treated as massless (MCFM) or massive (our calculation). The RHS plots show separately, 
for the massive case only, the scale dependence of the qq, gg and qg + qg contributions, as 
well as their sum. 
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Figure 3.10: Dependence on the renormalization/factorization scale of the rescaled difference 
between our NLO calculation (with mf, ^ 0) of the total Zhh cross section and the 
corresponding result computed using MCFM (with = 0) for the inclusive and exclusive 
cases (with /ir = /i/ = /io) respectively. The error bars indicate the statistical uncertainty of 
the Monte Carlo integration. 



is illustrated in Figures 3.11(a) and 3.11(b) for the inclusive and exclusive case, respectively. 
We see that the NLO QCD corrections affects the cross section quite substantially, in 
particular, for low values of the m^j invariant mass. In each figure the right hand side 
plot gives the ratio of the NLO and LO distributions. We stress the fact that the LO and 
NLO distributions are not just rescaled, which is clear from the RHS plots of Figures [3.111 



Figures 3.12(a) and 3.12(b) compare the NLO da/dmf,i distributions obtained from the 
massive and massless bottom-quark calculations. The results with = have been 
obtained using MCFM. As expected, most of the difference between the massless and massive 
bottom-quark cross sections is coming from the region of low m^j invariant mass, both for 
the inclusive and exclusive case, where the cross sections for mb ^ are consistently below 
the ones with m\, = 0. This is better emphasized in the right hand side plots, where we 
show the ratio of the two distributions, da{mb 7^ 0)/da{mh = 0). For completeness, we 
also show in Figure 13.131 the comparison between massive {mh 7^ 0) and massless {mh = 0) 
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Figure 3.11: The distribution da{pp Zhh)/d'm^,i in LO and NLO QCD. The right hand 
side plot shows the ratio of the LO and NLO distributions. 
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Figure 3.12: The inclusive and exclusive distributions da{pp Zhh)/dmi,i derived from our 
calculation (with ^ 0) and from MCFM (with = 0). The right hand side plot shows 
the ratio of the two distributions, da{mf, ^ 0)/da{mi) = 0). 
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Figure 3.13: The LO distribution da{pp Zbb)/dmi,i derived from our calculation (with 
rrih 7^ 0) and from MCFM (with rUb = 0). The right hand side plot shows the ratio of the 
two distributions, dcr{mb ^ Q)/da{mb = 0). 



calculations at LO in QCD. The LO m^^; distribution for massive bottom-quarks has been 
obtained both from our calculation and from MCFM, which implements the ^ option 
at tree level, and both results agree perfectly. In general, mass effects are similar at LO and 
NLO. To illustrate this in more detail we show in Figure [3. 141 the rescaled difference between 
the m^i; distributions obtained with our NLO calculation (with 7^ 0) and with MCFM 
(with rrib = 0) defined as in Eq. (13. 2p . We notice that, in the Zbb case, finite bottom-quark 
mass effects are relevant up to values of the m^^ invariant mass around 50 GeV. 

For ongoing searches of a light SM Higgs boson, regions with small m^j invariant mass 
are of relevance, as in many cases only one 6-jet is tagged (semi-inclusive studies) in order 
to increase the experimental statistics. In such studies m^^ invariant mass distributions are 
produced by using the leading two jets in the event, one of them being the only tagged 6-jet. 
In such case, a possible signal can come from a real emission in which the bb quark pair is 
recombined into a single b-jet (that is the partonic m^j invariant mass is small) and the extra 
light parton is seen as another jet. 
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Figure 3.14: The mf,5 distribution of the rescaled difference between our NLO calculation 
(with rrih 7^ 0) and MCFM (with = 0) for the inclusive (upper plot) and exclusive (lower 
plot) pp — ^ Zbb production. 
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CHAPTER 4 



CONCLUSION 



We have presented a full review of our calculation of NLO QCD corrections to W/Z bb 
production at hadron colliders including full bottom-quark mass effects [8^ [92] . We have 
shown results for total cross sections and bb invariant mass (m^j) distributions at the Tevatron 
Fermilab collider. We have found that, for such collider, the NLO QCD corrections reduce 
considerably the dependence on factorization and renormalization scales, in particular when 
considering exclusive cross sections where exactly two 6-quark jets are tagged in the final 
state. This then reduces the theoretical uncertainty of the cross section, from about 40% 
at LO to 20% and 10% at NLO for the inclusive and exclusive cases, respectively. Even 
more importantly, we have found that NLO corrections change considerably the shape of LO 
distributions; that is, the NLO distributions are not simply a rescaling of those at LO. 

We have systematically compared our results to a calculation that considers massless 
bottom quarks and have found that this approximation overestimates the total cross section 
by about 8% for Wbb production and 10% for Zbb production. The mass effects are 
particularly relevant in regions with small rrif^i invariant mass. On the other hand, the 
massless calculation shows very similar dependence on the factorization and renormalization 
scales, as including mass effects is expected to affect very mildly such dependence. 

4.1 Outlook 

As we have stressed, our results are of relevance to the search for a SM-like Higgs particle 
in the VH {V = W, Z) associated production channel and to the measurement of single-top 
production, both processes of great interest to the high energy physics community. The low 
invariant mass region, where bottom-quark mass effects are most relevant, is important 
when a light SM-like Higgs particle {Mh ~ 100 — 140 GeV) is searched semi-inclusively and 
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by tagging only one 6-quark, in order to increase the experimental statistics. The low m^g 
invariant mass region is of course relevant for single-top production since in this case the 
whole rrif^i spectrum is relevant, as the kinematics of the process is broader and semi-inclusive 
searches are essential. 

We are currently studying the impact of our calculation on searches for single-top 
production, where we also consider final states with fewer than two 6-quarks. We study 
modifications in the total cross sections and implications for 6-tagging efficiencjl^]. 

The next natural step is to implement our calculation for the LHC and study its 
phenomenological impact. Since at the LHC gluon initiated processes are enhanced, we 
expect some fundamental differences to appear. In particular, q{q)g initiated subprocesses 
will play a bigger role and, given their tree level nature, will increase the dependence 
on renormalization and factorization scales and, to some extent, increase the theoretical 
uncertainty of the cross sections. 

Finally, we want to emphasize that the calculation performed can be naturally extended 
to other important processes, such as '-fti production, which might be studied at the LHC 
and could give a direct measurement of the electric charge of the t-quark. In the same 
direction, we can study 766 production at NLO in QCD. This process has considerable 
phenomenological implications, as it can put direct constraints on the 6-quark PDF, which 
so far has only been derived from the gluon PDF evolution. 

We can also study the associated production of a pseudo-scalar with heavy quarks, which 
is of relevance to searches for physics beyond the SM, particularly models with an extended 
Higgs sector, like supersymmetric models. 

On the theoretical side, we expect to extend the checks performed on box coefficients 
by using generalized unitarity, to triangle, bubble and tadpole coefficients. Even more, we 
can extract analytically expressions for the contributing rational pieces. This will be of 
importance in the development of efficient new techniques for their extraction, as they will 
represent a playground to, for example, recursion relation techniques when massive external 
and internal particles are present in multi-leg processes. 



^ Since the completion of this Dissertation, we have shown in Ref. [98] exphcit results for NLO QCD 
corrections to Wb production both at the Tevatron and at the LHC. 
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APPENDIX A 



Standard Model of Particle Physics 



The Standard Model of particle physics is a quantum field theory based on the gauge groups 
SU3 for color, SU2 for weak isospin, and Ui for hypercharge, as dictated by the local 
gauge symmetry invariance observed in the behavior of fundamental particles. The color 
quantum number is associated with the dynamics of the strong interactions, which by itself 
is the subject of Quantum Chromodynamics (QCD), while the weak isospin and hypercharge 
quantum numbers are fundamental to the dynamics of electroweak interactions. 
The SM Lagrangian can be written as 

^SM = C,YM + Cf + Ch + -^^Yuk , (A.l) 

where Cym is the Yang-Mills Lagrangian, Cf the fermion Lagrangian, Ch the Higgs 
Lagrangian and /^Yuk contains the Yukawa interactions of the theory. Cym describes the 
dynamics of the gauge fields (kinetic terms + self-interactions) and includes the following 
terms 

Cym = Cqcd + Cj^ + £y 

8 3 

= -i E ^r-^'^"'^ - 4 E ^i^'^^'"^ - lB,.B^r (A.2) 

a=l i=l 

The color field strength tensor is given by 

Gl, = d,Al - d^Al + g^r'^A^Al , a, 6, c = 1, . . . , 8 , (A.3) 

with y4j^ the eight color gauge fields (so called gluons), gi the dimensionless strong coupling 
constant and f^^^ the structure constants of SU^. Analogously, the weak isospin, F^^, and 
hypercharge, -B^jy, field strength tensors are given by 

Flu = d,Wl - d^W; + g^e'^'W/^Wj; , k = 1,2,3, (A.4) 
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B^, = d^B, - d,B^ , (A.5) 

where PF^ and B^ are the 4 electroweak gauge bosons (a hnear combination of which will 
become the weak W^^ and Z° weak gauge bosons plus the photon A^, as shown in Eq. (1A.16P ). 
g2 is the dimensionless weak isospin coupling constant and e^^'' are the structure constants 
of SU2. 

Throughout the body of this dissertation we have denoted the strong coupling gi as Qs, 
and the weak isospin coupling (72 as gw We also use the conventional definition ag = g1/ (47r). 

The second part of the SM Lagrangian in Eq. flA.ip . Cf, describes the fermion fields and 
their interactions with the gauge bosons. The fermion fields are classified as quarks, which are 
triplets under the color gauge group, and leptons, which have no color. Taking into account 
the fact that the W boson couples only to left-handed helicity states of quarks and leptons, 
this part of the Lagrangian is built such that right-handed and left-handed components of the 
fermion fields couple independently to the gauge bosons. Using the notation {SU2,SU^)y, 
to denote weak isospin, color, and hypercharge quantum number assignments of the fermion 
fields, we can write that a quark weak doublet, Ql = is a (2,3)j^j and a quark weak 

singlet, u^, is a (1, Z)y^ (and similar for d/j). On the other side, a weak doublet of leptons, 
Li = (1^)^, is a (2, and a lepton weak singlet, e^j, is a (1, l)^^. The fermion Lagrangian 
Cf then can be written as 

+LLa^'V^LL + eRa^V^eR + -- - , (A.6) 

where the dots stand for similar terms for the remaining quarks and leptons. In Eq. (1A.6I) 
(7^ are the Pauli matrices {a^ = 1), and are the covariant derivatives corresponding to 
each field, 

V,Ql = {d, + gi'-A';^\'' + g2'^Wy + gs'-yiB,)QL, 
V^UR = (d^ + gi^A'^X + g3'^y2B^)uR , 
V^dR = {d^ + g,^A';^X^ + g-s'-ysB^)dR, 
V,Ll = id^ + g2'-Wy + g,'-y,B^)LL, 

'D^cr = {df, + g3'-y5B^)eR , (A.7) 
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where is the dimensionless hypercharge couphng constant, and r* and A" are the Pauh 
and Gell-Mann matrices for SU2 and SU^ respectively. 

Notice that a mass term for the fermion fields and for the vector boson fields (as needed 
for the weak vector bosons and is not allowed by gauge invariance. The last two 
terms in the Standard Model Lagrangian shown in Eq. flA.ip are introduced to remedy 
this problem. Indeed, the simplest way to preserve the gauge symmetry of the SM while 
generating massive electroweak gauge bosons is the so called Higgs mechanism, which we 
explain in the following. A separate step needs to be taken to introduce massive fermions, 
and we will discuss this below. 

The Higgs mechanism, in its simplest version HOOj llUll I1U2] , starts by adding to the 
model another field, called the Higgs field H, which transforms as a weak isospin doublet, a 
color singlet, and it has hypercharge yh- 

H = (^f^ ~ (2, 1),, . (A.8) 

Its dynamics is dictated by the Ch term in Eq. (lA.ip . which can be written as 

Ch = {V,H)\V^H) - V{H) , (A.9) 

where V^H = {d^ + 5'2|W^/^t' + gz^yhB^)H is the covariant derivative of H and V{H) is the 
most general renormalizable potential invariant under SU2 x f/i, 

V{H) = ix^H^H + \{H^Hf , (A.IO) 

with /i^ and A real parameters. A is a dimensionless parameter. 

If /i^ < the field configurations that minimize the potential V{H) has to satisfy: 

HLHvac = ^ = y • (A.ll) 

So, once /i^ < the Higgs field develops a vacuum expectation value (VEV), which is 
degenerate over the sphere defined in the last equation. Picking one configuration breaks 
this degeneracy, causing the vacuum of the theory not to be SU2 x Ui symmetric anymore. 
To illustrate the consequences, let us choose: 
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Indeed one can verify that this choice breaks the original gauge symmetry: 



SU2 xUi^ Ur , (A. 13) 

where Uf''^^ is the electromagnetic Ui symmetry. 

When Ch is expanded in the vicinity of the chosen minimum, by shifting the Higgs field 
as follows 

Ch becomes the Lagrangian of a real scalar field with mass rrih = 2f ^A, the physical Higgs 
boson. Moreover, a mass term for the gauge bosons is generated by the first term in Eq. flA.Qp . 
coming from 



(A.15) 



2" 2' 

The corresponding mass eigenstates, i.e. the physical gauge fields, are obtained by diagonal- 
izing the mass matrix of the vector fields Wj^ and 5^. The EW gauge bosons and Z^, 
as well as the photon A^, are expressed as: 

1 



-g,B, + g2Wl 



+ 9l 

92Ba + 93W. 



3 



A, = , (A.16) 

V ^2 + 9i 



with the associated masses: 



M2 
w = l92'v . 

Ml = \{9l + 9l)v\ 

Ml = . (A. 17) 

These simple relations are found to agree with experiment {My/ = 80.4 GeV, Mz = 
91.2 GeV) with v ^ 174 GeV. This is, of course, approximate as they are results based 
on the classical, or leading order level of the theory. 

Finally let us focus on the last part of the SM Lagrangian presented in Eq. (lA.ip . This 



term, Cyuk, couples massive fermion fields to the Higgs field via Yukawa type interaction. 
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For example, the gauge invariant Yukawa coupling of the Higgs boson to the down quark, 
d, is 

-XdQiHdji + h.c. , 

where Q^^ = (u, d)/, and is the Yukawa coupling for the down quark. After the shift of 
Eq. (1A.14P this term gives the effective coupling 

-Ad^Qlf ^^^dR + h.c, (A.18) 



'V2 \v + h 
which gives a mass term to the down quark with 

= ^ , (A.19) 
V 2 

and defines the coupling between the down quark and the physical Higgs particle to be 
-\d/V2. 

Similar terms are added for each massive fermion field. Then Cynk will contain 9 arbitrary 
parameters, the Yukawa couplings, standing for 6 quark masses and 3 lepton masses. This 
completes the classical SM Lagrangian, which is the main theory used throughout this 
dissertation. Several extensions of the SM follow the same prescription and ultimately 
break the EW symmetry spontaneously by introducing several Higgs fields. On the other 
hand, alternative mechanisms to break the EW symmetry have also been proposed like 
"Technicolor", where scalars are seen as strongly bound states of fermions, "Little Higgs" 
models, where more complex global symmetries of the fundamental fields are considered to 
justify the existence of a relatively light Higgs, or more speculative theories that explain 
the breaking of the EW gauge symmetry in terms of extra dimensions. However, most 
alternatives to the Higgs mechanism encounter phenomenological difficulties and till now 
the most successful theory in describing experimental data is still the SM with a weakly 
coupled Higgs boson. 

Feynman Rules 

From the Lagrangian one can readily extract the Feynman rules of the theory. Now we 
present the set of rules we have employed throughout our calculation. 

Let us start by writing the propagators for quarks and gluons, this last one in the Feynman 
gauge: 
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"^dfJ-iy cab 



(A.20) 



+ m) 



P 



p2 _ jj^2 _|_ 



{A.21] 



For the coupling between quarks and gluons and gluon self-interactions we have the 
vertices: 



a,iJ, 




(A.22) 



o 
o 
o 
o 



b,iy 




(A.23) 




-Wi [j J [g g - g g , 

_|_ jace f'^'^^ (^g^^^ gP'^ — g^" g^^ 
^fadejbce^g,ugpa_g,pgu. 



(A.24) 



where 2t"j are the Gell-Mann matrices and / the structure constants of SU{3). 

The weak boson couplings to fermions are given by the vertices (stripped of color indices): 



W,iJ. 




2V2 



7^^(1-75)1///. 



(A.25) 
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cosOw 



(A.26) 



where Vffi are the entries of the CKM mixing matrix, and the vector, (yfy, and axial, (7^, 
couphngs for the Zff vertex are given by: 

gl^ = ^Ti - sine jQf, (A.27) 

/ _ _^Tf 

9 A — 2 

where sin^vi' = 93/ ^ di, if the third component of the weak isospin matrix for the 
fermion /, and Qf is the charge of the fermion. 
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APPENDIX B 



Tree Level Amplitudes for W/Z bb Production 



In this Appendix we present explicit expressions for the tree level amplitudes that appear 
in the full calculation of W/Z bb production. This is a straightforward application of the 
Feynman rules shown in the last section of Appendix |Al As we mentioned there, we denote 
the strong coupling gi as Qs, and the weak isospin coupling g2 as gw 

Tree level amplitude for qq' — > Wbb 

The contributing tree level Feynman diagrams are shown in Figure 12.21 
Given the momenta assignment: 



where Vqq> is the CKM mixing matrix, e^(pvv-) is the polarization vector of the W boson and 
we denote by v and u the spinors for the external fermionic fields (for more details look at 
the Feynman rules in Appendix |A]) . 

Tree level amplitude for qq — > Zbb 

The contributing tree level Feynman diagrams are shown in Figures 12.21 (with V = Z) and 
12.101 for subprocess qq Zbb with the Z weak boson emitted from initial and final fermion 
lines respectively. 



the LO amplitude can be written as: 



Aoiqq' -> Wbb) 



= W. 




(B.l) 
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Given the momenta assignment: 



the LO amphtude can be written as: 

Ao{qq Zhh) = ig] ^'^ e* {j)z) , ^'^'^ 



{qi-PzY 

2 9w 



COS 9v 



1-^2 +Pz 

l^{9l + 9ll^)u, 
9up 



+ iib7 



i>b + i>z + mb 
[{Pb+PzY - ml 



nib 



[{-Pb-Pz)^ - ml 



Y{9t + 9''.l^)v-b 



iB.2) 



where gl and g{ are the vector and axial couphng constants for fermion /, e^{pz) is the 
polarization vector of the Z boson and we denote by v and u the spinors for the external 
fermionic fields (for more details look at the Feynman rules in Appendix lAl). 

Tree level amplitude for gg — ^ Zbb 

The amplitudes Aq^s, Ao,t, and Ao,u introduced in Section [2.4.11 can be written as: 

(^f^iQl) ^A(l2)(^*p{pz)ubA!^'^fv-b , 
(^Ml)(^,y{(l2)^*p{Pz)UbA'^'^/vi , 

eM(9i)e<.fe)e;(pz)Mfe^^Jufe , (B.3) 

where ^^{pz) is the polarization vector of the Z boson, v and u are the spinors for the external 
fermionic fields, and Aq^^ , Aq^, and Aq'I,^ represent the total s— channel, t— channel, and 
M— channel amplitudes, corresponding to the diagrams in Figure I2.11[ More explicitly: 

A 

A. , 



Ao,s 


= 


9w 


Ao,t 


= wl 


9w 


Ao,u 


= 


9w 



pup 
0,s 




A{2)pup 




pup 
Q,t 




A(2)pup 


A{'i)pUp 


pup 
0,u 


_ aWp-i^P 
~ •^0,u 


A{2)pup 


I A{'i)^^^p 
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where 

s (Pb + PzY-mi 



with 



mi] 



[{pb + Pzr - '"-bi 
s^'^iin + VzY-ml]^ 

\KPh - qiY - [(^2 - pir - ^bi 

^i'"^" = X?'^^''(/i-^,gi-g2) , 

At^'" = 51-^2) , (B.5) 



are the individual amphtudes for the s— channel, t— channel, and u— channel diagrams in 
Figure Em 
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APPENDIX C 
Scalar Integrals for W/Z bb Production 



In this Appendix we present a collection of all the one-loop IR-divergent scalar integrals 
that appear in the W/Z bb NLO QCD calculation. For completeness we will also include the 
UV-structure of the UV-divergent scalar integrals, namely tadpoles and self-energies. We 
will use the following notation for the scalar m-point integral: 



/0(gi 



4-d 



dH 



m. 



m—l. 



1 



{2Tif [t^ - m2][(t + - m?] . . .[{t + + ■ ■ ■ + qm-if 



m: 



id] 



which corresponds to the topology illustrated in Figure IC.ll We have denoted by t the loop 
momentum, by {qi} (i = 1, . . . ,m, qm = — Yl^=i 1j) incoming momenta to the 

diagram and by {rrii} (z = 0, . . . , m — 1) the set of masses corresponding to the propagators 
in the loop. The integration over the loop momentum is performed over d = 4 — 2e (e = 




Figure C.l: Topology of one-loop Feynman integrals. We denote by t the loop momentum, by 
{qi} the set of incoming momenta and by {rrii} the set of masses in the internal propagators. 
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unless stated differently) dimensions. We will denote integrals 10 with m = 1,2,3,4 and 5 
denominators by AO, BO, CO, DO and EO respectively, as conventional. 

Most of the expressions presented in this Appendix exist in the literature (see for 
example [103^ 17511771 11U41I5U] ). except for the box integrals labelled in the following as 11.2 and 
11.3, which as far as we know have not been presented explicitljjl. We have calculated these 
integrals using two independent techniques and found agreement. For IR-finite integrals we 
have used expressions from [103j and cross checked them with the FF-package [57] . 

The expressions presented below are correct up to 0{e), since this is what is needed 
for a full NLO calculation. We include only the real pieces □ of the integrals (according to 
the kinematics we use) as the imaginary parts do not contribute to the NLO partonic cross 
section. We will use the kinematics of the process: 

^(gi) + j{q2) v{py) + 6(pb) + h{pi) , 

with pv = Qi + <l2 — Pb — Pb ^^'^ the on-shell conditions Q'l = = and = = m^. 
Although in this dissertation we treat the V vector boson as on-shell (i.e. = M^), the 
integrals in this Appendix are correct also for for the case of an off-shell V boson, as we 
always eliminate pv by using conservation of momentum. Therefore they can be used also 
if one had to include leptonic decays of the V boson. Moreover, all the following results are 
expressed in terms of the invariants Sij = {qi + qjY, where qi and qj are two of the external 
momenta. 

Scalar integrals can be related by rotation or reflection of the integration momentum: 

m—l 

/0(gi, . . . , g^.i; mo, ... , m^_i) = /0(- ^ q^, qi, . . . , qm-2; m^.i, mo, . . . , mm-2) 

i=l 

= /0(gm_i, . . . ,gi,m„_i, . . . ,mo) . (C.2) 

In the following we therefore give the minimal set of scalar integrals appearing in our 
calculation. 

Finally we will always factor out of the integrals the factor: 



^b= { ^] r(l + e) . (C.3) 



ml 



^ Since the completion of this Dissertation a full set of IR-divergent scalar integrals have been presented 
in Ref. |105j . and we found agreement with the corresponding box expressions presented there. 
^By real pieces we mean the real pieces of the integrals after factorizing i/(167r^). 
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AO integral 

AO and BO integrals are the only scalar UV- divergent integrals. The 1-point integral is 
IR-finite, and it is given by: 

AOim,) = -i^A4 ml (— + 1) . (C.4) 

We notice that AO{m — > 0) = 0, as it can be inferred from the fact that such dimensionful 
integral cannot be built from any kinematic invariant. We mention that the only other 
1-point integral appearing in our calculation is AO{mt), which can be obtained from the 
previous expression replacing nib and Afb by nit and Aft respectively. 

BO integrals 

All BO integrals are UV-divergent and their corresponding UV-pole part is: 

SO(p;mo,mi) = — . (C.5) 

UV-pole Ib'K'^ euv 

The only IR-divergent BO integral is the one with zero internal masses and with a light- 
like {q\ = 0) external momentum BO{qi;0,0). Indeed, there is no invariant available to 
build this integral, so it should vanish. One can understand this vanishing as a cancellation 
between the UV and IR behavior of this integral of the form: 

BO{q,;0,0) = -^K ( — ~ —) , (C.6) 

where we have made explicitly the UV or IR nature of the two single poles. Notice that the 
overall constant in front of the integral is a matter of convention. 

IR-divergent CO integrals 

Let us organize the integrals by the number of internal masses. 

I. All internal masses equal to zero 

1. CO(gi, 0, 0, 0), two on-shell massless legs: 

C0(„. 0. 0. 0) = —A/-, _ ^- + - 1„ + - -±--y (C.7) 
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2. CO{qi, —ph — Pb'i 0, 0, 0), one on-shell massless leg: 



C0(gi,-p6-p5; 0,0,0) 
1 /I 



167r2 



2 qi ■ i-pb-Pb) 



P12 



In 



mt 



-Pb-PbY, 



P12 



TX 



-Pb-Pbf, 



{G.i 



notice that for this case, in our kinematics, Si2 < so we have included absolute values 
in the logarithm arguments to avoid negative arguments. 

II. One non-zero internal mass 

1. CO(gi, —pb\ 0, 0, rrib), One on-shell massless leg and one on-shell massive leg: 



CO(gi,-pfe;0,0,mfc) 



167r2 



-Mb 



Si2 - ml I 2e2 e 



In 



mt 



Si2 - mf 



+;inMi 



Sl2 



mt 



Lio 



Sl2 



+ ln 1 



^ In 1-^ 



TT 



(C.9) 



2. CO(gi, q2 — Pb'i 0, 0, rUb), One on-shell massless leg: 



C0(gi,g2 -P6;0,0,mfe) 



167r2 



tA4 



In- 



mt 



ei^'^si2-"^b 1(^2 -Pb)^ - "2^1 



+ ln^ 



Si2 - mf 



In^ 



In 



mt 



\{q2-PbY-ml 



In 1 



(g2 - Pb)^ - m^ 



Li2 1 



\{q2-PbY -ml\ 



Si2 - ml 



+Li2('"^^"^'r^^^' ]-V 



7 



mt 



(C.IO) 



III. Two non-zero internal masses 



1. CQ{—pb,Pb + Pb't 0? ^b), Two on-shell massive legs: 

i 1 



CO(-pb, pfe Pi; 0, mb, m;,) 



167r2 



Mb 



PiPb + PbY 



1, 



:/2 - tt' 



(C.li: 
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where we have defined: 



/3 



4 ml 



and 



TT^ + 21n/51n 



-2Li2 



2/3 y V 2/3 

IR-divergent DO integrals 

I. All internal masses equal to zero 

We write the massless box expressions following the results of [77]. 

1. il>0(gi, —Pb — Pb; 0, 0, 0, 0), two adjacent on-shell massless legs: 



D0(gi,g2,-P6-K;0,0,0,0) 



167r2 



1 

+ - 

-2Li2 ( 1 - 



+ 



"■§23 



mlsi2 



[Pb + Pbf{-qi -q2+Pb+ Pbf 
' -qi- q2+Pb + PbY 



S23 



- 2Lio 1 - 



ml 



'~qi -q2+Pb + PbY 



-qi - q2+Pb + p-bf 

Sl2 



-In' 



TT 

~6 



-S12 \ 
-S23) 

(C.12) 



2. DO{qi, —pb — Pbi Q2'i 0, 0, 0, 0), two opposite on-shell massless legs: 



D0(gi,-pb-p5,g2;0,0,0,0) 



167r2 



■Sl2'S23 



m, 



2 \ ^ 



-S12 



m 



2 \ e 



"■523 



mf 



-Ml 



-Ml 



-2 Li2 ( 1 ^ - 2 Li2 1 ^ - 2 Li2 1 ^ 1 - 2 Li2 I 1 

S12 J \ S23J V S12 



M' 



Ml 

S23 



Ml Ml, , 
+2 Li2 ( 1 1 - In^ 

■5l2'523 



-■Sl2 
"■523 



(C.13) 



with M| = {pb + p-^y and M| = (gi + q2 - Pb - Plf- 
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II. One non-zero internal mass 



1. Z)0(gi, "Pbi 0, 0, 0, TTib), two adjacent on-shell massless legs and one one-shell 
massive leg: 



D0{qi,q2,-Pb,0,0,0,mb) 
2 In 



3 

2^ 



mt 



tA4 

.2 



1 



167r2 S12(S23-"^^) 



mt - S23 



+ ln!!^-ln 



mt 



si2 {qi + q2 - PiY - ml _ 



-2 In 



mt 



ml - S23 



mt 



In ( ^ 1 - In" 



mt 



{qi + q2- PiY - ml 



-2Li2 1 + 



(gi + q2-Pby 



mt 



ml - S23 



(C.14) 



2. Z>0(gi, — gi + Pb + Pb? —Pb't 0, 0, 0, m;,), one on-shell massless leg and two one-shell 
massive legs: 



^0(gi, -qi +Pb + Pi, -Pi\ 0, 0, 0, m^) 



tA4 



1 



167r2 Si2(s23-mg) 



1 1 

2?^7 



In 



In^ 



-I- In 



mt 



-Sl2 



mt 



-Sl2 



In^ 



+ In 



mt 



ml - S23 



-In 



mt 



"1^ - 523 



-M| + gi2 
"^b - 523 



In 



[m 



2^2 



- 2 In ( 1 



'12 



-_SV2_ 

-Ml 



In 



m^ - S23 



In 



m^ - S23 



A _ -M|(-5i2 + mg- 523 + M|) \ 
/ V - S23)(-Si2) / 



In I \ 1^ - m,^ + S23 



ml - S23 



-Lio 



-Li^ 



Lis 1 



In 



ml - 523 



-M|(-3i2+m^-S23 + M| 
(m^ - S23)(-Sl2) 



-M| + Si2 



m^ - S23 



-M| + gi2 
ml - S23 



Li5 



In 1 



- g23 
-M| + Si2 

-A/| - mg g23 

-Sl2 



Li2 



- g23 
-M| + Si2 



Li2 



"^g - g23 



(C.15) 



where we have defined the mass square of the second leg as M| = {—qi + Ph + PiY 
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3. DO{qi, — gi — ^2 + Pb + Phi Q2 — Phi 0? 0? 0, rn.5)|f|, one on-shell massless leg and one 
one-shell massive leg: 



DO{qi, -qi -q2+Pb + Pi, 12 - Pi', 0, 0, 0, m;,) 



1 

+ ln2 



In 



mt 



-S12 



+ In 



mt 



ml - S23 



In 



167r2 



mt 



1 



Si2[S23 - ml) 



mt 



S23 



mt 



4 



ml ) 3 



+2 In ( In 1^ (-^12 + M|)(m^-. 23 + M| J ^ _ 2 / " 



m 



S12IS23 - ?Tif 



mt 



In 



S23 + M. 



S12 



"^h / \ mt- S23 



+2U2 ( 1 
-2Li2 ( 1 



;-Si2 + M|)K-S23 + M| 



S121S23 - "ife 



-mlMl-{s23-ml){Mi-ml) 

2Li, ( 1 - !^i^^^i±M 
S12 



-Li2 1 



-Si2 + M| 

'Tifc - S23 

(ml - S23)M| 



g23 - ml 
Ml 



mlMl - (523 - ml){Ml - mf ' ' 

Ml + si2)M. 



Uo 1 



mlMi + {s2s-ml){Ml-ml) 



-U2 1 



^1 



mli-Ml + 512) - (S23 - ml){Mi - ml) 



+ Li2 1 



m 



2 -Ml 



mt 



+ 



Li2 (1 



-S12 + M|)M| 



V mli-Ml + S12) - (S23 - ml){Mi - ml) 
where we have defined M| = {—qi — q2 + Pb + PiY ^'^d M| = (g2 — Ps)^- 



(C.16) 



During the completion of this manuscript a compilation of one-loop scalar integrals has appeared on 
http://qcdloop.fnal.gov/, written and maintained by R.K. Ellis. The box integrals presented in sections II. 2 
and II. 3 of this Appendix are presented in there with one expression, corresponding to i?g in their notation. 
We have compared analytically the pole structure and numerically the finite real part of our expressions II. 2 
and II. 3 with the expression of -Bg, and wc found perfect agreement. (Sec footnote [1] in this Appendix). 
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III. Two non-zero internal masses 



1. DO{qi, —pi,q2;0,0,mh,mh), two opposite on-sliell massless legs and one on-shell 
massive leg: 



DO{qi, -pi, q2; 0, 0, nib, nib) 



1 



167r2 isi2-ml){s23-ml) 



1 1 



In 



ml - S23 



In 



mt 



In 



mt 



ml-si2j \{qi + q2 - PbY - ml 



2 S23 



mt 



In^ 



2 mb- Si2 



mt 



-21n(^£Hl±(?l±^i^)ln 



- Si2 



ml - si2 



-21n(^£i^±i?l±^^^lln 



ml - S23 



ml - S23 - Si2 + (gi + q2- VlY 

ml - S23 

ml - S23 - si2 + (gi + g2 - Pb^ 



-2U2 
-2U2 

-2U2 



ml - S23 - .§12 + (91 + g2 - Pb)^ 

- 512 

ml - g23 - ^12 + (91 + g2 - Pb)^ 

ml - S23 

{-S23 + (gi + g2 -Pb)^)(-gi2 + (gi + q2-PbY) 

{ml - S23){ml - S12) 



(C.17) 



2. Z)0(gi, — P5, — gi+Pb+PB; 0, 0, rrib, 7715), one on-shell massless legs and two opposite 
on-sliell massive legs: 



DO(gi, -pi, -gi +pb + Pb\ 0, 0, mfe, mb) 



1 



1 1 

^ + - 



In 



mt 



ml - S12 



In 



mt 



ml - S23 



mt 



ml - S12 



+ ln^ 



mt 



2 i^l - Si2\ 27r2 

+ m 



ml - S23 



3 +^"H^)+^"Hr 



m,l - S23 
(C.18 



with 



^^"2! {-qi+Pb + Pi)^ 
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3. D0{q2, —Pbi —Qi — Q2 -\- Pb -\- Pb'i 0? 0, m^, 7715), one on-shell massless leg and one 
on-shell massive leg: 



D0{q2, -pi, -qi -q2+Pb + Pi', 0, 0, nib, nib) 



1 



167r2 {si2 - ml){s23 - ml 



1 1 



In 



ml - S12 



In 



mt 



S23 - ml 



In 



mt 



-M2 



mt 



ln2 f "^fc ~ \ _ [ml - M'i 



mt 



+2 ln(^^i^)lnf^^^^^l+2 1n 



ml - S12 



ml - M2 



■§23 - mt 



mt 



-Si2 + Ml\ ^ ml - Ml 
ml - S12 



-2 Li2 — „^ ^^^^ — ^ - 2 Li2 



"^b - S12 



S23 - ml 
S23 - Ml -ml + S12 



S23 - ml 



.2Li2 ( K-^l)(^23-M|-mg + Si2) . _ 

{S2z-ml){ml- S12) ° 



where M| = (gi — pfc)^ is the mass square of the fourth leg, 



(C.19) 



In 



mt 



-M2 



\ In / ^3 



m. 



mt 



+ ln 



- Li2 



mt - S12 



mt 



Ml 



In 



-Lig 



[ml - Ml) - A+(-5i2 + Ml) 



ml - S12 



X4-SU + MI 



m. 



-M2 



-S12 + Ml 



mt 



Ml 



\+{-Si2 + Ml) + ml - Ml 



+ (A+ ^ A_ 



(C.20) 



Mi 



—Qi — C[2 + Pb + Pb) is the mass square of the third leg and 



A. 



i± Wi 



Ami 
Mi 



(C.21) 
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IV. Three non-zero internal masses 



1. D0{—ph,q2,—q2 + + Pb\^i''^bii^bii^b), One on-shell massless leg and two 
adjacent on-shell massive legs: 



DO{-ph, q2, -q2 +Pb + Pb, 0, rrih, m^, rrib) 



1 



1 fl-p 



167r2""" S23{si2 - ml)P 
2 In (Xs) In {ml/ {ml - S12)) + 2 In (Xs) In (l - X^) 
-7rV2 + Li2(A',A',)+ln2(A'3) 

-2 (Li2(A',A'3) + In (X,) In (1 - A'.A'g) + In In (1 - ) 
-2(Li2(A',/A'3) + In (Xs) In (1 - XJX^,) - In (A'3) In (1 - XJX^,) 



(C.22) 

where, using M| = (— ^2 +P6 + PiY the mass square of the third leg, we have defined: 



-3 = t/l-^ 

S23 



P3 



4 ml 



and 



1-/3 
1 + P 

1-/33 



2. DO{—pb, —qi — q2+Pb+Pbi 91 + ^2; 0, mt, m^, 7715), Two adjacent on-shell massive 
legs: 



DO{-pb, -qi -q2 + Pb + Ph Qi + <l2\ 0, m?,, m^,, mf,) 



tA4 



1, (^-13 
- In 



where we have defined: 



167r2 (si2-m2)s23/3 U \l + (3 
/9 = 



, (C.23) 



4 m? 



S23 
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Defining M| = (-gi -q2+Pb + Pi) , ^3 = (Qi + ^2)' 



P2 = X I 



4 ml 



2 



= 

X2 = 



4 



^3 

1+13 ' 
1-/^2 

1 + /32 ' 

and 

;^ _ 1-/53 

we can express the finite piece Xq as follows: 



2 In (A',) In (l - X^) - 2 In (Xs) In ^/(m^ - S12)) 
+7rV2 + LialA-^) + ^^(A'a) + In^X^) 

-U^iXgX^X^) + Li2(^.A'2/A'3) + U2{XsX:,IX2) + Li2(A',/A'2/A'3) 

+ In (1 - X.X^X^) In (A",) + In (1 - A", A's A's) In {X^) 
+ In (1 - ;f, A'2A'3) In (A'3) + In (1 - A'. A'2/A'3) In (A'.) 
+ In (1 - A',A'2/A'3) In {X2) + In (1 - A',A'2/A'3) In (l/A'3) 
+ In (1 - A", A'3/A'2) In (A",) + In (1 - A", A'3/A'2) In (l/A'2) 
+ In (1 - A',A'3/A'2) In (A'3) + In (1 - A',/(A'2A'3)) In (A",) 

+ In (1 - A',/(A'2A'3)) In (l/A'2) + In (1 - A',/(A'2A'3)) In (l/A'3) 

(C.24) 

Reduction of box and pentagon scalar integrals 

An interesting property of m-point scalar integrals in d = 4 — 2e dimensions is that they can 
be expressed as a linear combination of (m — l)-point scalar integrals in 0? = 4 — 2e plus the 
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corresponding m-point scalar integral calculated in d = 6 — 2e dimensions. This has been 
proved in Ref. [75] and is summarized in the following formula: 



+ (m - 5 + 2£)c„4«-='' 



i=l 



(C.25) 



where the integral on the left hand side corresponds to (— l)'""^^/0(gi, . . . , qm-i', ^o, ■ ■ ■ , fnm-i), 
in the notation of Eq. fIC.ip . Jm"^"^'^'' denotes the same m-point integral calculated in 
d = 6 — 2e dimension and /^-i correspond to the (m — l)-point integrals obtained from 
Ijn by taking out the i-th denominator and are calculated in ci = 4 — 2e dimensions. 



The Cj 



coefficients in Eq. (1C.25I) are given by: 



Co 



i=l 



(C.26) 



where the matrix S',-,- is built from invariants as: 



1 ( 2 

2 ^ * 



4) • 



(C.27) 



with pI = (g^ + . . . + qj)^. 

We have used Eq. (1C.25P twice. First, using Eq. (1C.25P we have cross checked our 
calculation of the IR-divergent box integrals II. 2 and II. 3 that did not exist in the literature. 
Second, we have used it to compute the pentagon scalar integrals. In the case of i?0-functions. 



Eq. flC.25p is particularly useful because the coefficients in front on /, 



(d=6-2e) 



is of (9(e) and 



therefore does not contribute (given that /^'^-6-2^) is IR-finite). According to Eq. (^C25\\ . 
i50-functions are therefore calculated as the linear combination of five scalar integrals in 
= 4 - 2e. 
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APPENDIX D 



Reducing Tensor Feynman Integrals 



At one-loop, Feynman diagrams with m external legs may contain Feynman tensor integrals 
with up to m denominators of the form: 



Jn^^-^"(gi, . . .,qm-i,mo, ■ ■ ■ .rrim-i) = 



(27r)'' [t^ - ml][{t + - mj] ■ ■ ■ [{t + + ■ ■ ■ + qm-if 



'IJ " " " LV' "T yi "T ■ ■ ■ -r Lirn-lj — 'm-m-li 

which corresponds to the topology illustrated in Figure ID.li In Eq. fID.ip t is the loop 
momentum, c? = 4 — 2e is the dimension of the loop momentum integration, Hi, . . . , fi^i are 
Lorentz indices associated to each power of loop momentum appearing in the numerator, 
{li} {'i = ^, ■ ■ ■ i"^, Qm = — Yl^=i ^j) is if incoming external momenta connected to 

the loop and {rrii} (z = 0, . . . , m — 1) is the set of masses associated with each propagator 
in the loop. The index "n" is associated with the rank of the tensor integral. As a short 



(D.i: 




Figure D.l: Topology of one- loop Feynman Integrals. We denote by t the loop momentum, 
{qi} the set of incoming momenta and {rrii} the set of masses in the internal propagators. 
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hand notation, depending on the number of denominators, or external momenta, we will 
write I = A for one leg, I = B ioi two legs, and so on. Integrals with no powers of loop 
momentum in the numerator are called "scalar" (see Appendix [U]) and are labelled by n = 0. 
So, for example, a scalar triangle integral will be denoted by CO, while a rank 4 box integral 
by DA. 

One notices immediately that the tensor integral In^^'"^^ is a symmetric tensor and it 
can be expressed as a linear combination of (symmetric) tensor Lorentz structures built from 
the external momenta {gj} and the metric tensor, g^'^ , if the set {qi} is not complete. In this 
case (i.e. for m > 5) g^'^ can be expressed as a linear combination of momenta. 

The coefficients of the Lorentz structures are rational functions of scalar integrals and 
invariants built from the external momenta {g^}, the masses {rrii} and the integration 
dimension d. For example, a straightforward computation of the m = 2, n = 1 tensor 
integral gives: 



where we have defined the coefficient of the only Lorentz structure as B^^\qi,m,m) = 
— |BO(gi, m, m). 

The systematic reduction of tensor integrals to a linear combination of Lorentz structures 
is known as the Passarino-Veltman (PV) method [72], and consequently we will call the tensor 
integral coefficients Jn-PV functions, depending on the topology and rank of the integral. 
In the following we will write explicitly the tensor integrals we have encountered in our 
calculations and give an example of how they are recursively reduced (for a review see for 
example [73]). 

Starting with 2-point integrals, we write (suppressing the arguments of the PV functions): 



i?l^(gi, m, m) = — i?0(gi,m, m 



(D.2) 



Bl^{qi,mQ,mi) 
B2^'{qi,mo,mi) 



(D.3) 
(D.4) 
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For 3-point integrals we write: 

Cl(gi,g2,mo,mi,m2r = C^gf + ^(2)^^^ , (D.5) 



C3(gi,g2,mo,mi,m2)^'^^ = C^^^^^g^^q^^ + perm) + C^^^^\g^"'q^ + perm) 

+C^^^^\qUi<l2 + perm) + C^^^^\qi^q^2<il + Perm) , (D.8) 

where the term "perm" accounts for all terms, obtained from permutations of Lorentz indices, 
that are needed to obtain a symmetric tensor. 
Finally for 4-point integrals we write: 

Dl(gi,g2,g3,rno,mi,m2,m3)'^ = D^^^gf + D (2)^2^ + D (2)^3^ , (D.9) 
D2(gi, 52, gs, mo, mi, m2, ms)'^^ = D^^^^g^' + Z^^^^g^r + D^^''\W2 + ^^''^ga^gg + 

D(i2)(grg2^ + grg2") + D^'''\qUl + g^ga") + (D-10) 
^^''Hg2%^ + g2g3") , 

^3(gi,g2,g3,mo,mi,m2,m3)'^"'' = L'^^^'H^^'^V + Perm) + + perm) + 

/)(003)(^Mi^^p + perm) + (D.U) 

i^^^^^gfgrg? + ^^^^^^^g^^g^g^^ + ^^^^^^gs^gsgf + 

D'^^^^\qUi<l2 + perm) + D'^^^'^\q^q'(ql + perm) + 
D'-^^^Hq2(l2(li + perm) + D^^^^\q^q^q(; + perm) + 
D^^^^\q^qlq{ + perm) + D'^^^^\q^qlqP^ + perm) + 

^^'''Hgfgs gf + perm) , 
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DA{qi, ga, gs, rriQ, mi, ms, m^yp" = 

/^(oooo) ^^M.^p<x ^ pg^^^ _^ ^ ^(oon) ^gi^-qPq- + perm) + 

je{i,2,3} 

J2 D^''^^\g^''{qPqJ + qPq^)+ perm) + + 

jje{i,2,3} i6{l,2,3} 

E ^^^^^^■^(grgr9fgJ + Perm)+ D^^^^^^q^q-^q^q^ + perm) + 

i,je{i,2,3} ije{i,2,3} 

J2 D^''''\qUU';ql + perm) . (D.12) 



i,i,fc6{l,2,3} 



Naturally, the calculation of pentagon diagrams (like the ones shown in Figures 12.61 12.151 
and I2.19P involves up to £'4-PV functions, but, as explained in Sections 12.2.21 12.3.11 and 
12.4.21 we have only used directly these functions for cross checks, due to the large numerical 
instabilities associated with them. Their structure can be easily inferred from previous PV 
functions, by just omitting terms with g^-^ tensors, and will not be presented here. 

Expressions for the Jn-PV functions can be obtained recursively by saturating the tensor 
Lorentz indices of Eq. (ID. 10 and Eqs. (]D.3P - (lD.12p with the same external momentum or 
gP-'^ tensor. Using relations like 

t-qi = \[{{t + q^Y - ml) - {t" - ml) - (qf + ml - mj)] , (D.13) 

allows one to simplify the loop momentum dependence in the numerator of Eq. (ID.ip against 
some of the denominators, thereby relating higher rank m-point tensor integrals to lower rank 
and/or less external point ones. 

Comparing the coefficients of the reduced tensor structures obtained by saturating 
Eq. fID.ip and the corresponding expressions in Eqs. flD.3p - flD.12p with the same external 
momentum or g^''^ tensor provides a set of equations that fully determines the coefficients of 
the tensor integrals represented in Eqs. flD.3p - flD.12p . 

Let us consider as an example the case of DA^'^p'^ . The set of D4-PV functions is one 
of the most demanding irreducible pieces that appear in the W/Z bb NLO QCD calculation 
(without taking into account the reducible E-PV functions). They are needed to account 
for the abelian contributions of the box diagrams 

b[Y'^^ to Af, shown in Table EH 
To write the set of equations that allow the expression of all 22 D4-PV functions in terms 
of lower rank functions, let us denote by Xp'^'-* the term proportional to the tensor structure 
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QjQk^i from the contraction of DA^'^^" in Eq. fID.ip with qi^ {i,j,k,l = 1,2,3), that is for 
example: 

X(^i2) _ coefficient gfgj'g^ of in DA^'^P^^q^^ . (D.14) 

We will denote by xf'^'^^ the coefficients of the structure g^'^q'l obtained from the contraction 
of DA^"'P'' with g^^. 

More explicitly, here is what one obtains by contracting DA^^'"^ of Eq. (ID.ip with the 
external momenta (i = 1, 2, 3): 

D4^"''"'{qi,q2,q3, mo, mi, 1712,1713) qi„ = ^ 

C3'"'^(gi + q2, qs, mo, ma, mg) - C3'"'''{q2, qs, rrii, mg, mg) 

+q'^ C2'"P{q2, qs, mi, m2, 1713) + C2''^(g2, gs, mi, 1712, m^) + gf C2'"'{q2, gs, mi, m2, ms) 
-gf gi Cl^(g2, ga, mi, m2, mg) - gf gf Cl''(g2, ga, mi, m2, mg) 
-gfgi Cl^(g2,g3,mi,m2,m3) 

+Q1Q1Q1 C0{q2, ga, mi, m2, mg) - (gi ■ gi + mg - mi) D3^''''(gi, g2, gs, mo, mi, m2, mg) 

(D.15) 

DAf"'P''{qi,q2,q3,mo,mi,m2,m3) q2a = ^ 

C3^''^(gi, g2 + gs, mo, mi, ms) - C3''''^(gi + g2, gs, mo, m2, ms) 

-(92 ■ q2 + 2g2 ■ gi + m? - ml) DS'^'^iqi , g2,gs,mo,mi,m2,ms) . (D.16) 

Di^^P^i^qi, g2, gs, mo, mi, m2, ms) gsa = ^ 

C3''''''(gi, g2, mo, mi, m2) - C3^''''(gi, g2 + gs, mo, mi, ms) 

-(gs ■ gs + 2gs ■ (gi + g2) +m2 - mg) D3^''^(gi, g2, gs, mo, mi, m2, ms) , (D.17) 

where all the m-point tensor integrals used are defined in Eq. fID.ip . Expressing the Dn and 
Cn tensor integrals in Eqs. (1D.15P - (1D.17P using Eqs. (lD.5p - (]D.lip . one reduces them to linear 
combinations of symmetric tensor structures with three Lorentz, built of external momenta 
and the metric tensor, whose coefficients are precisely the X^-'''^^ defined before Eq. f lD.Mp . 

Eqs. f lD.i5p - flD.i7p in this form have them to be compared to the expressions obtained 
from Eq. flD.12p when contracting DA^'^p" in there with gj^- {i = 1,2,3). The result of the 
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comparison can be cast into the following matrix formula: 

qi ■Q2 ql g2 ■ gs X 
qi ■ gs g2 ■ gs ql j 

£){0011) ^(0012) ^(0013) ^(1111) /)(2221) ^(3331) ^(1112) ^(1113) ^(1122) ^(2213) ^(1133) 
2^(0012) £)(0022) £)(0023) £)(1112) ^(2222) £){3332) £){1122) £)(1122) £)(2221) £){2223) £){3312) 
£)(0013) £)(0023) £)(0033) £)(1113) £){2223) £){3333) £){1123) £)(1123) £)(2213) £){2233) £)(3331) 

/ j(001)_^(0000) j(OOl) j(OOl) j(lll)_3£){0011) 

^(002) j(002)_^(0000) j(002) j(222) 

^ j(003) j^003) j(003)_^(0000) j(333) 

j(lll) j{112)_2£)(0012) j(113) _ 2^(0013) 

j(222)_ 3^(0022) j{222) j(112) _ ^(0011) j{113) 

-^^333) x(333)_ 2^(0033) jai2) jai3) _£)(0011) 



X(221) _£){0022) j{223) j{331) _^(0033) \ 

^(221) _2^(0012) j(223)_2^{0023) j{331) 

j{221) j^223)_^(0022) j^331) _2£)(0013) ^ 



(D.18) 

where X^'^''* is defined before Eq. (lD.14p . The solution of this equation determines almost 
all 22 D4-PV coefficients in terms of lower-rank, lower-point PV functions. More accurately, 
Eq. (ID.lSp can be solved for all the coefficients Z}(*-?'^'), except for D^^'^^'^\ In fact, all other 
coefficients will depend on it. To solve completely all Z)4-PV functions, we need then another 
relation for Z^C^ooo)^ This relation can be obtained by contracting DA^''^'^ with gp„ and 
applying a relation like Eq. fID.lSp . The result is: 

DA^''P''{qi, q2,q3,mo,mi,m2, m^) = 

C2^''{q2, qs, mi, 1712, m^) - q^ CV{q2, qs, ^i, m2, m-s) - q'[ Cl^{q2, qs, mi, 1712, m^) 

+gfgr C0(g2,g3,"2i, 7712,7713) 

+ml D2'"'{qi,q2,q3,mo,mi,m2,ms) . (D.19) 

Notice that all terms on the RHS of the previous equation depend on C77i-PV functions 
(77i < 4), and D2-PV functions. From this relation one them obtains the desired relation for 
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^(0000)^ namely: 

D^''''^ + E 2* ■ + (rf + 2) = gs, mi, m2, mg) + ^(0°) , 

(D.20) 

after which all the 22 D4-PV functions are determined in terms of lower rank and lower point 
PV functions, kinematic invariants built from the external momenta and internal masses as 
well as the dimension of the momentum space d. 

When the same procedure is repeated for all the PV functions, one ends up with a full 
reduction to scalar integrals and rational functions depending on the external momenta, the 
internal masses and the dimension d. 

Finally, it is evident from Eq. (ID.ISP how the Gram determinants (GDs) that we discussed 
extensively in Section I2.4.2[ appear in the PV-reductions. The GD for 4-point functions 
is indeed the determinant of the first matrix on the LHS of Eq. (ID.lSp . The solution 
of Eq. flD.181) involves at least one inverse power of it for each coefficient D^'^^^^\ Even 
two powers can appear in those cases when some entries of the RHS depend on Z}(0'^*-?) 
(z, j = 0, 1, 2, 3). So at least one power of GD appears in the reduction of D4-PV functions 
to D3-PV functions and other lower rank and lower point coefficients. By the time they are 
completely reduced to scalar integrals, the D^'^^^^^ tensor integral coefficients will contain at 
least four power of the 4-point GD. 
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APPENDIX E 



Using Quadruple Unitarity Cuts to Check Coefficients 

of Scalar Box Integrals 



In this Appendix we review a set of non-trivial cross checks performed on pieces of our 
calculation of NLO QCD corrections to W/Z bb production including full 6-quark mass 
effects. By using generalized unitarity cuts, specifically quadruple cuts as presented by 
Britto, Cachazo and Feng (BCF) [23] , we have been able to check results for the coefficients 
of scalar box integrals, which we have computed analytically using Passarino-Veltman (PV) 
reduction (see Appendix [D|) . 

The Calculation 

We show now explicit results for the quadruple cut involving a top loop, which contributes to 
the virtual corrections to the subprocess gg — > Zbb (see Section I2.4.3P shown in Figure IE.lt 
with the kinematics: 

9{qi)9{q2) ^ Z{p,)b{pb)b{pi), ql = ql = pi - ml = pi -ml=pl-Ml = 0. (E.l) 

This quadruple cut has the feature that only the B^} Feynman diagrams with a top loop 
(see Figure I2.18P contribute to it, that is all pentagon diagrams in the subprocess vanish 
under the cut. We notice that this box topology is of particular interest because it involves 
the most intricate irreducible tensor integrals of the computation, namely D4-PV functions, 
as described in Appendix [Dl 

To extract coefficient of the scalar box b[^^ , from the quadruple, we simply use a 
similar approach to the example given in Section 3 of Ref. [23]. We get, analogously to 
BCF's Eq. (3.1), the expression: 

= E - ^t) {{i + qif - ml) ((£ + + q,f - ml) ((£ + p,f - ml) < [ , (E.2) 
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Figure E.l: Topology of the example presented. It corresponds to two Feynman Diagrams 
given by the two possible orientations of the fermion line. The loop in the fermion is a top 
quark (with mass rrit). 



where i± corresponds to the two solutions of the on-shell conditions: 

{i\f = ml {i + q,f = ml {i + q, + = ml {£ + = m^,} , (E.3) 

which we solve (following BCF [23]) by using the parametrization: 

i = aq^+ (3q2 + ap, + pP^ , = e^'"'"' qi.q^pPz.- (E.4) 

The term B^] in the RHS of Eq. (lE.2p corresponds to the expression of the Feynman 
diagrams where the loop integral is frozen, i.e. the one obtained by substituting the loop 
momentum with the solutions I = l±. We show the solutions to the on-shell conditions in 
Eq. (lE.Sp in the following subsection. 

We have compared analytically both results for the box coefficient S^]- , from our standard 
computation and from the generalized unitarity technique, and they agree. We do not write 
such expressions, as they are quite cumbersome and not too illuminating. 

We have performed similar checks for other coefficients of scalar boxes. In general, we 
have calculated, using PV reduction, all diagrams in such way that all tensor integrals are 
reduced to a set of 1-, 2-, 3- and 4-point scalar integrals. So we can actually extract tadpole, 
bubble, triangle and box coefficients, as well as rational functions, from any given set of 
Feynman diagrams. They can be used as a playground for on-shell one-loop techniques like, 
for example, the ones discussed in Ref. [2]. 
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Solutions to on-shell conditions 



Here we write explicitly the solutions to the on-shell conditions in Eq. (lE.Sp . by using the 
parametrization shown in Eq. flE.4p . One obtains a set of four quadratic equations for a, /3, 
(T and p which can be solved straightforwardly. The two solutions l± are given by: 



a 



/3 



a 



^ [vl^ + 2Vi'iV22. + 2t;i4t;23 + 2^24 ^^23 - ^'23^^12 - ^'23^^34 " ^23"^^ + 2i;i3W24 



-^^24^^34 - Vi2ml + + 2Wi4i;24 - ^24"^^ " f 12^'34 " ^24^^12 ) / ^ , 
1 (^^24 + ^^23 - "^6 - ^^34)(^^12 - ^^13 " ^m) 



±- 



2 
1 

--Vi2 

2 

'2 A 



(f 24 + ^^23 - ml- f34) 

A 



:e.5) 



for the two possible signs in the p expression. We have defined: 



and 



{vi2vi-i + 2V2-imlvi'i - 2t;23t;i2f34 + 2V2Zmlvu - 2f23^'l2"^^ + 2t;23f24t'l2 
+2m^W24^^13 - 1^12^34 - 2^24^12^^ - 2mfvi2ml + t;i2^^24 - 2l'24^'l2W34 



+2vi2mlv34 + vuml + VuV^^^ + 2m^VuV24 ) / Vu 



A = P4 ■ -P4 = ^^13^^23 + ^^13^^24 + f 14^^23 " V^ml - V12V34 + t;i4f 24 



(E.6) 



(E.7) 



which corresponds to the Gram determinant of the process. The invariants Vij are defined 
as Vij = qi ■ qj, with qi and q2 as in Eq. (lE.ip . ^3 = pi, and q^ = pi- 

We mention that the box coefficient in Eq. flE.2p depends only on even powers of p, giving 
then a solution which is a rational function of invariants as expected. 
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APPENDIX F 



Phase Space Integrals for the Emission of a Soft Gluon 
in the two Cutoff PSS Method 



Phase space soft integrals for qc^ initiated W/Z bb production 

In this appendix we collect the integrals which we have used in calculating the results in 
Eq. (12.361) starting from Eq. (12.341) . For a more exhaustive treatment of the formalism used 
we refer to Refs. [791 EH], from which the results in this appendix have been taken. 
We parameterize the soft gluon (i-momentum in the qq' rest frame as: 

= . . . , sin6'i sin6'2, sin6'i cos6'2, cos6'i) , (F-1) 

such that the phase space of the soft gluon in (i = 4 — 2e dimensions can be written as: 



io 



(F.2) 



Then, all the integrals we need are of the form: 



Jo Jo 



{a + bcos9i) ^ 



(yl + S cos 6^1 + C sin 6^1 cos 6'2)' 



(F.3) 

In particular we need the following four cases. When ^ B"^ + C'^, and b = —a, we use 
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(dropping terms of order O {{d — 4)"^)): 



while when 6 7^ —a we use: 

/(o,i) 



TT 



a{A + B) {d~A 



+ In 



(A + B) 



In^ 



A + B 



- -In^ 
2 

+ 2 Lis 

- 2Li2 



A + VB^ + 
A-y/WT&^ 

B + VB^ + C2 



A + B 



TT 



In 



Li2 



(d-A) 



A + + C2 

2V52 + C2 



A + V52 + C2 



+ -In 



4 V A - VB^ + C2 



(F.4) 



(F.5) 



/(0,2) 



27r 



A2 - 52 - C2 



X 



2^ ^ V52 + C2 I A - V52 + C2 



Finally, when ^2 = 52 _|_ (^2^ ^^^^ ^ _ have: 



27r 
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Phase space soft integrals for gg Zbb 

In this Appendix we collect the phase space integrals for a final state soft gluon that are 
used in calculating the results reported in Eqs. fl2.8ip and fl2.84l) . We parameterize the soft 
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gluon d-momentum in the gg rest frame as shown in Eq. (IF.ip . We have seen that the phase 
space of the soft gluon in (i = 4 — 2e dimensions can be written as in Eq. (lF.2p . 
Then all the integrals we need are the following four: 
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(F.8) 



where we have used the set of kinematic invariants in Eq. f l2.23p . (3bi and A^g are defined in 
Eq. ([226]), K and A^^ after Eq. ^IM^ and Mb in Eq . fl239|) . Moreover we have denoted by 
F(j)i,pf) the function: 



F{Pi,Pf) = In^ 
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where is the angle between partons i and / in the center-of-mass frame of the initial state 
partons, and 

/ o 

Sif 



Pf 



'I 



(3 f cos 9if 



(F.IO) 



All the quantities in Eq. ( IF. 91) can be expressed in terms of kinematical invariants, once we 
use Sif = 2pi-pf and: 



pI 



s-s-bv + ml 



2yi 



and 
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2yi 



(F.ll) 



with SfY = {pf +PvY {V = Z). 
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